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Dear  Readers: 

This  issue  marks  the  third  year  of  publication  for  the  Building  Research  Journal.  It  represents 
the  technology  transfer  mission  of  the  Building  Research  Council  of  the  University  of  Illinois  at  Ur- 
bana-Champaign  and  the  National  Consortium  of  the  Housing  Research  Centers. 

This  journal  has  gone  through  a  remarkable  maturing  and  refinement  process  over  the  last  four  is- 
sues. During  this  time,  we  have  broadened  the  scope  to  include  Infrastructure  Planning  and  Develop- 
ment topics.  With  the  Fall  1993  issue,  we  started  a  new  feature  section  titled:  "News,  Updates,  and 
Seminars."  This  section  was  originally  designed  to  provide  the  opportunity  to  share  and  discuss 
ideas  and  events  in  an  informal  format.  With  the  current  issue,  we  have  expanded  the  "News,  Up- 
dates, and  Seminars"  section  to  include  selected  activities  of  the  members  of  the  National  Consor- 
tium of  the  Housing  Research  Centers. 

In  past  issues,  we  have  tried  to  balance  the  papers  so  that  they  represent  several  emphasis  areas 
of  the  journal.  In  this  issue,  we  have  deviated  from  the  practice  and  decided  to  focus  on  one  area.  The 
four  papers  in  this  issue  represent  the  area  of  "Building  Design,  Operation,  and  Performance."  We 
would  like  to  receive  your  comments  in  this  regard. 

The  paper  by  G.  K.  Yuill  and  E.  D.  Werling  deals  with  the  calculation  of  residential  heating  en- 
ergy requirements.  The  analysis  of  one  of  the  main  causes  of  indoor  air  quality  problems  and  its 
spread  through  the  HVAC  system  is  presented  in  the  paper  by  Ossama  A.  Abdou  and  Francis  A. 
Sando.  Ardeshir  Mahdavi  and  Khee  Poh  Lam  discuss  building  performance  simulation  in  their  pa- 
per. The  final  paper  by  Elizabeth  Kossecka,  Jan  Kosney,  and  Jeffrey  E.  Christian  presents  a  mathe- 
matical model  for  estimating  thermal  resistance  of  a  wall. 

We  would  like  to  thank  the  contributors,  the  reviewers,  the  Building  Research  Council  staff  and 
the  members  of  the  Editorial  Board  for  their  input.  We  would  also  like  to  thank  the  readers  for  mak- 
ing this  publication  a  success.  Your  submissions,  subscriptions,  comments  and  suggestions  are  wel- 
come. 
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Interpolating  Residential  Heating  Energy  Requirements 

G.  K.  Yuill  and  E.  D.  Werling 


ABSTRACT 


INTRODUCTION 


Life  cycle  cost  optimization  of  residential 
heating  equipment  requires  knowledge  of  an- 
nual heating  energy  requirements.  Tradition- 
ally, the  degree-day  method  has  been  used. 
More  accurate  methods,  using  computer  simu- 
lation, are  available,  but  they  require  hourly 
weather  data  which  is  not  available  at  many  lo- 
cations. This  problem  could  be  solved  by  inter- 
polation between  the  nearest  sites  with 
available  hourly  weather  data. 

BLAST  was  used  to  calculate  heating  energy 
requirements  for  three  houses  at  34  sites  in  the 
northeastern  United  States  with  available 
hourly  weather  data.  Then,  24  of  the  34  sites 
were  used  to  evaluate  the  accuracy  of  various 
prediction  methods,  including  the  degree-day 
method,  distance  weighting  factors,  and  linear 
and  non-linear  interpolation.  Ten  boundary 
sites  were  needed  for  interpolation  at  sites 
near  the  boundaries  of  the  region  considered. 

The  accuracy  of  the  degree-day  method  was 
worse  than  any  of  the  other  methods.  Standard 
errors  of  four  to  seven  percent  were  observed 
for  the  better  interpolation  methods,  compared 
with  over  11  percent  for  the  degree-day 
method.  The  best  interpolation  methods  allow 
the  analyst  to  predict  residential  heating  en- 
ergy requirements  at  a  site  for  which  weather 
data  are  not  available,  without  introducing  sig- 
nificant error. 


Dr.  G.K.  Yuill  is  a  Professor  of  Architectural  Engineering  at 
The  Pennsylvania  State  University. 

E.D.  Werling  is  completing  a  Master  of  Science  in  Architec- 
tural Engineering  at  The  Pennsylvania  State  University. 


Life  cycle  cost  optimization  of  residential 
heating  equipment  requires  a  knowledge  of  an- 
nual heating  energy  requirements  for  the 
house  of  interest,  local  utility  rates,  and  the 
costs  of  equipment  components  being  consid- 
ered. Utility  rates  and  equipment  costs  are  not 
difficult  to  determine.  However,  annual  heat- 
ing energy  requirements,  which  depend  on 
house  characteristics,  occupant  behavior,  and 
local  weather,  can  be  more  difficult. 

Traditionally,  the  degree-day  method  has 
been  used  to  estimate  heating  energy  require- 
ments. However,  it  is  not  accurate,  and  is  be- 
coming less  accurate  as  houses  become  more 
energy  conserving  and  appliance  energy  use  in- 
creases. More  accurate  methods,  using  com- 
puter simulation  of  building  energy  use,  are 
available,  but  they  require  hourly  weather 
data,  which  is  not  available  at  many  locations. 
This  problem  could  be  solved  by  interpolation, 
using  simulated  heating  energy  requirements 
at  the  nearest  sites  for  which  hourly  weather 
data  are  available.  This  paper  presents  several 
such  methods  for  estimating  annual  heating  en- 
ergy requirements,  and  compares  the  accuracy 
of  each. 

GENERAL  APPROACH 

In  this  project,  the  BLAST  computer  pro- 
gram (Building  Loads  Analysis  and  System 
Thermodynamics,  U.S.  Army  CERLXBLAST 
1991)  was  used  to  calculate  the  annual  heating 
energy  requirements  for  three  house  designs  at 
34  sites  in  the  northeastern  United  States  for 
which  the  required  weather  data  were  avail- 
able (Figure  l)(Table  1).  Then,  24  of  the  34  sites 
(all  but  the  boundary  sites)  were  selected  as 


10 


Building  Research  Journal 


Table  1 :  BLAST  Reference  Locations 


City 


Elevation 


Latitude 


Longitude 


Hartford,  CT 

19ft 

41.7 

Wilmington,  DE 

74ft 

39.7 

Washington,  D.C. 

14ft 

38.9 

Caribou,  ME* 

624ft 

46.9 

Portland,  ME 

43ft 

43.7 

Baltimore,  MD 

148ft 

39.2 

Patuxent  River,  MD* 

50ft 

38.3 

Boston,  MA 

15ft 

42.4 

Concord,  NH 

342ft 

43.2 

Lakehurst,  NJ* 

20ft 

39.9 

Newark,  NJ 

7ft 

40.7 

Albany,  NY 

275ft 

42.8 

Binghamton,  NY 

1590ft 

42.2 

Buffalo,  NY* 

705ft 

42.9 

Massena/Richards,  NY* 

207ft 

44.9 

New  York,  NY 

13ft 

40.7 

Rochester,  NY 

547ft 

43.1 

Syracuse,  NY 

410ft 

43.1 

Akron,  OH 

1208ft 

40.9 

Cincinnati,  OH* 

758ft 

39.2 

Cleveland,  OH 

777ft 

41.4 

Columbus,  OH 

812ft 

40.0 

Dayton,  OH* 

1002ft 

39.9 

Toledo,  OH* 

669ft 

41.6 

Youngstown,  OH 

1178ft 

41.3 

Allentown,  PA 

387ft 

40.7 

Erie,  PA 

731ft 

42.1 

Harrisburg,  PA 

308ft 

40.2 

Philadelphia,  PA 

5ft 

39.9 

Pittsburgh,  PA 

1137ft 

40.5 

Scran  ton,  PA 

930ft 

41.3 

Providence,  RI* 

51ft 

41.7 

Burlington,  VT 

332ft 

44.5 

Charleston,  WV* 

939ft 

38.4 

*Boundary  site 

source:  ASHRAE  Fundamentals, 

1993 

72.7 
75.6 
77.0 
68.0 
70.2 
76.7 
76.2 
71.0 
71.5 
74.1 
74.2 
73.8 
76.0 
78.7 
74.9 
73.8 
77.7 
76.1 
81.4 
84.5 
81.9 
82.9 
84.2 
83.8 
80.7 
75.4 
80.2 
76.8 
75.3 
80.0 
75.7 
71.4 
73.2 
81.6 
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Massena 


Hartford    \  ^\    a 

_^^i^Prov  idence 


Figure  1 .  Weather  Data  Locations 


test  sites  to  evaluate  the  accuracy  of  the  vari- 
ous prediction  methods.  The  accuracy  of  each 
method  was  quantified  by  comparing  the  pre- 
diction method  estimates  with  the  BLAST 
simulations  at  the  same  24  sites.  The  ten 
boundary  sites  were  needed  to  interpolate  heat- 
ing energy  requirements  at  sites  near  the 
boundaries,  but  they  themselves  could  not  be 
used  for  interpolation,  since  they  were  at  the 
edges  of  the  region  considered.  For  example, 
simulations  of  Ohio  and  West  Virginia  sites 
were  needed  to  interpolate  heating  energy  re- 
quirements in  western  Pennsylvania  locations. 

The  compared  prediction  methods  included: 

•     the  degree-day  method; 


weighting  computed  estimates  at 
the  nearest  three  sites  by  inverse 
distances; 

weighting  computed  estimates  at 
the  nearest  three  sites  by  inverse 
square  distances; 

linear  interpolation; 

second-order  interpolation; 

inverse  square  distance  weighting 
with  a  linear  correction  for  eleva- 
tion; 

linear  interpolation  with  a  linear 
correction  for  elevation;  and 


second-order  interpolation  with  a 
linear  correction  for  elevation. 
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BLAST  simulations  were  used  as  bench- 
marks rather  than  using  actual  measured  heat- 
ing energy  requirements  because  of  the  large 
number  of  locations  required  to  evaluate  the 
prediction  methods,  and  because  BLAST  has 
been  validated  extensively  (Andersson  et  al. 
1980)(Bauman  et  al.  1983)(Herron  et  al. 
1980)(Yuill  1985  and  1986)(Yuill  and  Phillips 
1981).  BLAST  gives  accurate  results  when 
building  and  weather  input  data  are  accurate. 

THE  BLAST  SIMULATION  MODEL 

General  Description 

Three  typical  modular  house  sizes  were  mod- 
eled using  BLAST:  a  1-story  1,500  ft2  (139  m2) 
house  [House  1],  a  2-story  2,250  ft2  (209  m2) 
house  [House  2],  and  a  2-story  3,000  ft2  (279 
m2)  house  [House  3].  The  simulation  models 
were  based  on  the  plans  for  a  modular  house 
with  2,250  ft2  (209  m2)  of  finished  floor  area  (ex- 
cluding basement)  built  over  a  full  basement 
foundation  in  1992  in  central  Pennsylvania 
(the  Base  House).  The  house  has  standard 
wood  frame  construction,  with  2x6  insulated  ex- 
terior walls,  insulated  first  floor,  insulated  sec- 
ond story  ceiling,  double  pane  low-emissivity 
(low-e)  windows,  and  a  vented  un-insulated  at- 
tic. A  summary  of  the  component  areas  and  the 
overall  thermal  resistances  is  given  in  Table  2. 

Each  house  was  modeled  in  BLAST  as  three 
thermal  zones  (the  living  space,  an  uncondi- 
tioned basement,  and  an  attic)  with  typical 
schedules  for  occupancy,  lighting,  and  electric 
appliances  for  their  respective  house  size.  Total 


annual  energy  consumption  for  lighting  and  ap- 
pliances (excluding  HVAC)  was  7,000  kWH  in 
House  2,  based  on  data  from  a  West  Penn 
Power  study  of  residential  electric  demand  in 
houses  of  the  same  size  (Burt,  Hill,  Kosar,  Rit- 
telmann  Associates  1992).  For  Houses  1  and  3, 
lighting  energy  was  scaled  to  floor  area  and  ap- 
pliance energy  use  was  scaled  to  occupancy. 
BLAST  simulations  were  run  with  TMY  (typi- 
cal meteorological  year)  or  TRY  (test  reference 
year)  weather  data  for  each  of  the  34  locations 
in  the  northeastern  United  States. 

Basement  Model  and  Ground  Temperatures 

No  available  hourly  energy  analysis  pro- 
grams have  accurate  algorithms  for  estimation 
of  basement  heat  loss,  since  it  would  not  be 
practical  to  simulate  three  dimensional  tran- 
sient heat  transfer  on  an  hourly  basis.  Instead, 
BLAST  calculates  basement  heat  loss  by  as- 
suming one-dimensional  conductive  heat  trans- 
fer for  the  basement  floor  and  the  below-grade 
portion  of  basement  walls.  It  also  assumes  a 
ground  temperature  for  each  month.  The  diffi- 
culty is  to  determine  an  effective  ground  tem- 
perature for  each  month  that  gives 
approximately  the  same  results  as  a  more  pre- 
cise, but  prohibitively  complex,  three  dimen- 
sional transient  heat  transfer  simulation. 

Mitalas'  procedure  for  calculation  of  residen- 
tial below  grade  heat  loss  (Mitalas  1987)  was 
used  to  determine  these  effective  monthly 
ground  temperatures  for  each  of  the  34  loca- 
tions in  the  BLAST  simulation  model.  Mitalas' 
procedure  uses  Basement  Heat  Loss  Factors 
for  calculation  of  basement  heat  loss,  based  on 


Table  2.  House  Component  Areas  and  R-Values 

House  1 

House  2 

House  3 

Component 

Total  Area 

Total  Area 

Total  Area 

R-Value* 

Walls 

1067  ft2  (99  m2) 

1847  ft2  (172  m2  ) 

2215  ft2  (206  m2) 

19  (3.3) 

1st -Floor  Floor 

1500  ft2  (139  m2) 

1375  ft2  (128  m2) 

1500  ft2  (139  m2  ) 

19  (3.3) 

lst-Floor  Ceiling 

1500  ft2  (139  m2) 

30  (5.3) 

2nd-Floor  Ceiling 

1375  ft2  (128  m2) 

1500  ft2  (139  m2) 

30  (5.3) 

Windows 

164  ft2  (15.2  m2) 

264  ft2  (24.5  m2) 

328  ft2  (30.5  m2) 

3.125(0.55) 

2  Wood  Doors 

41ft2  (3.8  m2) 

41  ft2  (3.8  m2  ) 

41  ft2  (3.8  m2  ) 

7.19(1.27) 

Glass  Patio  Door 

40  ft2  (3.7  m2  ) 

40  ft2  (3.7  m2  ) 

40  ft2  (3.7  m2) 

3.45  (0.61) 

*  hr-ft2-°F/Btu  (m2- 

°C/W) 
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Table  3.  BLAST  Infiltration  Coefficients  from  LBL  Correlated  Daily  Infiltration  Rates 


CITY 


D 


Scranton,  PA 

1.3255 

-0.0128 

0.000408 

4.76  E-8 

Harrisburg,  PA 

1.3355 

-0.0127 

0.000407 

5.82  E-8 

Washington,  D.C. 

1.3541 

-0.0122 

0.000281 

1.39  E-7 

Portland,  ME 

1.3709 

-0.0129 

0.000324 

9.97  E-8 

foundation  type  and  geometry,  basement  insu- 
lation levels,  local  weather,  indoor  tempera- 
ture, and  soil  thermal  conductivity. 

Infiltration 

Infiltration  is  modeled  in  BLAST  with  a  cor- 
relation of  the  form: 

I  =  ImaxF[A  +  B(Ti-T0)  +  Cv+Dv2]    (1) 

where, 

/  =  estimated  infiltration  rate,  cfm  (Us); 

Imax  =  user-specified  maximum  infiltration  rate, 
cfm  (IVs); 

F  =  user-specified  fractional  infiltration  sched- 
ule, dimensionless  (1  for  this  study); 

Tj  =  average  indoor  temperature  for  the  hour 
(thermostat  setting),  °F  (°C); 

T0  =  average  ambient  outdoor  temperature  for 
the  hour,  °F  (°C); 

v  =  average  wind  speed  for  the  hour,  ftVmin  (m/s); 

A-D=  user-specified  coefficients. 

While  this  model  offers  user  flexibility  in 
specifying  infiltration  rates,  it  is  not  based  on 
physical  reality.  A  more  accurate  infiltration 
model,  based  on  the  Lawrence  Berkeley  Labora- 
tory (LBL)  model  (ASHRAE  1993)(Sherman 
1986)  takes  the  form: 


I  =  L^A(Ti-T0)+Bv' 


(2) 


where, 

L=  measured  Effective  Leakage  Area  (ELA), 

2  2 

in.   (cm  ); 

v  =  average  wind  speed  for  the  hour,  mph  (m/s); 

A  =  stack  coefficient,  cfm2/in.4oF  ((L/s)2/cm4oC) 
[0.0156  (0.000145)  for  House  1,  and  0.0313 
(0.000290)  for  Houses  2  &  3]; 

B  =  wind  coefficient,  cfm  /in.  mph 
{(L/s)2/cm4(m/s)2} 

[0.0039  (0.000104)  for  House  1,  and  0.0051 
(0.000137)  for  Houses  2  &  3]. 

To  accurately  model  infiltration  for  this 
study  the  ELA  of  the  Base  House  was  meas- 
ured using  a  blower  door  (CGSB  1986),  and  the 
average  infiltration  rate  was  calculated  for 
each  day  of  the  year  using  the  LBL  model. 
Weather  data  for  Scranton,  PA,  which  most 
closely  matched  the  climate  at  the  site  of  the 
house  tested  in  central  Pennsylvania,  were 
used  in  the  calculations.  The  daily  averaged  in- 
filtration rates  were  then  correlated  to  an  equa- 
tion of  the  same  form  as  the  BLAST  infiltration 
model,  to  solve  for  the  BLAST  infiltration 
model  coefficients,  A-D  (Table  3,  Stranton).  The 
correlation  was  excellent,  especially  at  high  in- 
filtration rates,  indicating  that  the  BLAST  infil- 
tration model,  with  these  correlation 
coefficients,  accurately  calculated  LBL  model 
infiltration  rates. 

To  check  the  validity  of  the  Scranton  coeffi- 
cients at  other  sites,  the  procedure  was  re- 
peated for  three  additional  locations.  The  new 
coefficients  (Table  3)  were  used  to  estimate  infil- 
tration in  the  Scranton  house  and  compared 
with  the  original  set.  The  correlations  for  all 
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three  cases  were  excellent  (standard  errors  of 
less  than  1  cfm),  showing  that  the  coefficients 
developed  using  Scranton  weather  could  be  ap- 
plied at  any  location. 

BLAST  Simulations 

Energy  analysis  results  from  the  BLAST 
simulations  are  shown  in  Table  4.  Annual  heat- 
ing energy  requirements  varied  from 
14,600,000  BTU  (4,280  kWH)  for  House  1  in 
Patuxent  River,  Maryland  to  67,000,000  BTU 
(19,600  kWH)  for  House  3  in  Caribou,  Maine. 
Winter  design  heat  loss  was  calculated  by  defin- 
ing a  BLAST  temporary  design  day,  using  the 
ASHRAE  design  conditions  (ASHRAE  1989) 
and  ignoring  solar  and  internal  gains.  Design 
heat  loss  ranged  from  13,500  Btu/Hr  (3.96  kW) 
for  House  1  in  Washington,  D.C.  to  33,600 
Btu/Hr  (9.85  kW)  for  House  3  in  Caribou, 
Maine. 

PREDICTION  METHODS 

Degree-day  Method 

The  degree-day  method  is  a  simplified  heat- 
ing energy  consumption  estimation  procedure 
(ASHRAE  1989)  based  on  the  assumption  that 
heating  energy  consumption  varies  linearly 
with  the  difference  between  the  mean  daily  out- 
door temperature  and  a  base  temperature. 
This  base  temperature  is  the  temperature  at 
which  solar  and  internal  gains  to  the  house  ex- 
actly offset  the  overall  heat  loss  through  the 
house  envelope  by  conduction,  infiltration  and 
ventilation,  and  is  assumed  to  be  65°F.  Consid- 
ering seasonal  equipment  efficiency,  the  degree 
day  method  can  be  used  to  estimate  heating  en- 
ergy consumption.  By  ignoring  equipment  effi- 
ciency, the  degree-day  equation  used  to 
estimate  the  annual  heating  energy  require- 
ment is: 


E  =  C, 


24qLDD 
AT 


(3) 


where, 


E  =  estimated  annual  heating  energy  require- 
ment, BTU; 

qL  =  design  heat  loss,  including  infiltration, 
BTU/Hr; 


DD  =  number  of  65°F  degree-days  per  year; 
A  T  =  winter  design  temperature  difference 

W  indoor  ~  •*  design  )>    *  > 

Cq  =  empirical  correction  factor  for  heating  ef- 
fect versus  65°F  degree-days. 

The  degree-day  method  is  easy  to  use  and  de- 
gree-days data  are  available  for  many  locations 
in  the  U.S.  and  Canada.  However,  the  method 
is  based  on  the  design  heat  loss  for  the  house 
which  does  not  directly  account  for  solar  and  in- 
ternal gains,  or  transient  heat  transfer,  and  re- 
lies on  a  simplified  estimate  of  infiltration 
losses.  These  simplifying  assumptions  can  re- 
sult in  large  errors  (Claridge  et  al.  1987)(Fis- 
cher  et  al.  1982). 

Variations  in  building  materials,  thermostat 
settings,  and  appliance  energy  consumption  af- 
fect the  base  temperature,  and  thus  the  accu- 
racy of  the  degree-day  method.  The  use  of 
variable  base  degree-days  (ASHRAE  1993)  to 
determine  the  actual  base  temperature  (usu- 
ally below  65°F)  for  a  house  would  improve  the 
accuracy  of  the  degree-day  method.  However, 
this  would  require  the  use  of  more  detailed 
weather  data.  With  detailed  weather  data  avail- 
able, more  accurate  simulation  methods  would 
be  preferred.  The  question  remains,  which 
method  to  use  when  detailed  weather  data  are 
not  available? 

Inverse  Distance  Weighting  Factor 

Weighting  factor  methods  were  used  to  esti- 
mate the  annual  heating  energy  requirement 
for  a  house  in  a  given  location,  by  using  simu- 
lated annual  heating  energy  requirements  for 
the  same  house  in  three  locations  where  hourly 
weather  data  were  available.  The  estimate  was 
made  by  weighting  the  "known"  heating  energy 
requirements  by  the  inverse  of  the  distance  be- 
tween the  chosen  location  and  the  reference  lo- 
cations. The  reference  locations  chosen  were 
the  nearest  sites  with  available  weather  data, 
so  that  they  formed  a  triangle  around  the  loca- 
tion of  interest  (Figure  2).  This  method  gives 
nearer  locations  greater  weight,  and  results  in 
a  straight  numerical  average  if  all  three  loca- 
tions are  equidistant  from  the  location  of  inter- 
est. 
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Table  4 

.  BLAST  Simul 

aled  Annual  Heating  Energy  Requirements  and  Winter  De 

sign  Loads 

CITY 

House  1 

Annual       Design 

(MMBtu)    (kBtu/H) 

House  2 
Annual       Design 
(MMBtu)    (kBtu/H) 

House  3 
Annual        Design 
(MMBtu)    (kBtu/H) 

Hartford,  CT 

25.4 

14.6 

35.7 

21.0 

42.7 

25.5 

Wilmington,  DE 

19.0 

14.3 

27.1 

20.9 

32.4 

25.6 

Washington,  D.C. 

17.1 

13.5 

24.7 

19.8 

29.3 

24.3 

Caribou,  ME 

40.4 

19.7 

55.9 

27.7 

67.0 

33.6 

Portland,  ME 

30.4 

16.5 

42.5 

23.6 

50.8 

28.6 

Baltimore,  MD 

18.4 

14.3 

26.5 

20.8 

31.7 

25.4 

Patuxent  River,  MD 

14.6 

14.2 

18.4 

19.6 

21.9 

24.1 

Boston,  MA 

25.1 

16.5 

35.8 

24.0 

43.2 

29.6 

Concord,  NH 

29.7 

16.9 

41.3 

24.1 

49.3 

29.3 

Lakehurst,  NJ 

18.2 

13.9 

26.0 

20.3 

31.0 

24.8 

Newark,  NJ 

19.2 

14.3 

27.3 

20.9 

32.8 

25.6 

Albany,  NY 

26.7 

16.8 

37.3 

23.9 

44.6 

29.1 

Binghamton,  NY 

34.1 

16.7 

47.5 

24.0 

56.8 

29.2 

Buffalo,  NY 

32.2 

15.7 

45.4 

22.7 

54.7 

27.7 

Massena/Richards,  NY 

35.3 

18.6 

48.7 

26.4 

58.4 

32.1 

New  York,  NY 

20.6 

14.8 

29.5 

21.7 

35.4 

26.7 

Rochester,  NY 

32.9 

15.7 

46.3 

22.8 

55.7 

27.9 

Syracuse,  NY 

31.5 

16.3 

44.3 

23.2 

53.1 

28.3 

Akron,  OH 

26.6 

15.5 

37.5 

22.4 

44.9 

27.4 

Cincinnati,  OH 

18.7 

16.1 

26.5 

23.3 

31.6 

28.6 

Cleveland,  OH 

28.5 

16.2 

40.0 

23.4 

47.8 

28.6 

Columbus,  OH 

22.6 

15.4 

32.0 

22.1 

38.3 

26.9 

Dayton,  OH 

23.2 

16.6 

32.8 

24.0 

39.3 

29.3 

Toledo,  OH 

27.6 

16.1 

38.8 

23.1 

46.5 

28.1 

Youngstown,  OH 

29.8 

16.3 

42.0 

23.4 

50.3 

28.6 

Allentown,  PA 

24.2 

15.5 

34.1 

22.4 

41.0 

27.4 

Erie,  PA 

29.5 

15.1 

41.7 

21.9 

50.1 

26.8 

Harrisburg,  PA 

20.1 

15.2 

28.5 

22.1 

34.0 

27.1 

Philadelphia,  PA 

18.9 

14.4 

27.0 

21.1 

32.3 

25.9 

Pittsburgh,  PA 

25.8 

15.8 

36.3 

22.7 

43.4 

27.7 

Scranton,  PA 

25.8 

15.5 

36.2 

22.3 

43.2 

27.2 

Providence,  RI 

24.7 

14.9 

34.9 

21.7 

41.8 

26.6 

Burlington,  VT 

32.6 

17.7 

45.3 

24.9 

54.2 

30.2 

Charleston,  WV 

16.6 

14.2 

24.0 

20.6 

28.4 

25.1 

min      14.6 

13.5 

18.4 

19.6 

21.9 

24.1 

max     40.4 

19.7 

55.9 

27.7 

67.0 

33.6 
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Location  2 

Qi 


Figure  2.  Weighting  Factor  Methods 

The  inverse  distance  weighting  factor  proce- 
dure used  to  estimate  heating  requirements 
can  be  represented  by  the  following  equation: 


Qi     v      l 


(4) 


where, 

Qi  =  BLAST  simulated  heating  energy  re- 
quirement for  a  reference  location,  BTU; 

Z?i  =  distance  between  the  location  of  interest 
and  a  reference  location,  degrees. 

Since  latitude  and  longitude  were  readily 
available  for  every  location,  the  distances  be- 
tween the  location  of  interest  and  each  of  the 
reference  locations  were  calculated  in  degrees 
from  the  following  equation: 


A  =  V(L-Li)2  +  (0-0;): 


(5) 


where, 


Lj  =  longitude,  degrees,  subscript  indicating  a 
reference  location; 

4>  i  =  latitude,  degrees. 

This  distance  formula  is  not  exactly  correct 
since  latitude  and  longitude  are  only  equiva- 
lent at  the  equator.  However,  for  ease  of  calcu- 
lation, these  errors  were  ignored,  since  all 
locations  were  subject  to  approximately  the 
same  error.  These  errors  did  not  affect  the  heat- 
ing requirement  estimates. 

The  inverse  distance  weighting  factor 
method  was  an  improvement  over  the  degree- 
day  method  because  it  was  based  on  detailed 
computer  simulations  of  the  house  in  question, 


which  accurately  considered  the  effects  of  solar 
and  internal  gains,  heat  storage  (thermal  tran- 
sients), and  infdtration  losses.  However,  the 
method  did  not  account  for  differences  in  cli- 
mate not  explainable  by  latitude  and  longitude 
differences,  such  as  the  effects  of  elevation, 
proximity  to  large  bodies  of  water,  etc. 

Inverse  Square  Distance  Weighting  Factor 

The  inverse  square  distance  method  used 
was  the  same  basic  procedure  as  the  inverse 
distance  method.  The  difference  was  the 
method  placed  greater  weight  on  the  nearest 
reference  locations  by  squaring  the  distances  in 
the  weighting  factor.  The  estimates  were  calcu- 
lated with  the  following  equation: 


n 

1    ui         l    ui 


(6) 


This  method  offered  a  slight  improvement 
over  the  inverse  distance  method  because  it 
gave  extra  weight  to  the  nearer  locations. 

Linear  Interpolation 

This  method  was  similar  to  the  weighting 
factor  technique  since  the  heating  energy  re- 
quirement estimate  was  based  on  known  heat- 
ing energy  requirement  values  from  nearby 
locations.  However,  instead  of  weighting  by  dis- 
tance alone,  a  linear  relationship  was  assumed 
between  geographical  location  and  annual  heat- 
ing energy  requirements.  That  is,  heating  en- 
ergy requirement  was  assumed  to  be  directly 
related  to  latitude  and  longitude.  Based  on  this 
assumption,  an  equation  for  heating  energy  re- 
quirement was  developed  using  the  known  val- 
ues from  nearby  reference  locations,  from 
which  the  house  annual  heating  energy  require- 
ment in  the  location  of  interest  was  estimated. 
The  equation  was  of  the  form: 


E  =  ah  +  b<&  +  c 


(7) 


where  the  coefficients  a,  b,  and  c  were  deter- 
mined by  solving  the  equation  simultaneously 
with  the  data  points  from  the  reference  loca- 
tions. With  three  data  points  (locations  with 
known  annual  heating  energy  requirements) 
an  exact  solution  to  this  system  of  three  linear 
equations,  was  found  using  algebraic  methods. 
When  more  data  points  were  used,  it  became  a 
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simple  multiple  linear  regression  and  was 
solved  by  the  method  of  least  squares. 

Second-Order  Interpolation 

This  method  was  based  on  the  same  princi- 
ple as  the  linear  interpolation  method,  except 
that  a  higher  order  relationship  (non-linear) 
was  assumed  between  location  and  annual 
heating  energy  requirements.  The  relationship 
took  the  form: 


2  2 

E  =  ah  +  60  +  aZJ0>  +  dL  +eO  +/ 


(8) 


This  method  required  data  from  six  or  more 
locations  to  solve  for  the  six  coefficients. 

Number  of  Weather  Data  Locations 

The  number  of  weather  data  locations  used 
in  the  linear  and  second  order  interpolation 
methods  was  varied  to  determine  its  effect  on 
the  accuracy  of  the  interpolations.  It  was  ex- 
pected that  increasing  the  number  of  locations 
would  improve  the  statistical  accuracy  of  each 
method.  The  minimum  sample  size  that  gives  a 
good  prediction  would  be  desirable  for  most  cir- 
cumstances. The  base  case  for  each  interpola- 
tion method  was  the  minimum  number  of 
locations  required  to  get  a  solution  (i.e.  three 
for  linear  and  six  for  2nd  order  interpolation). 
The  largest  number  of  locations  used  was  24, 
the  number  of  locations  with  verifiable  BLAST 
simulations.  In  these  large  sample  correlations, 
three  separate  correlations  were  performed  (24 
sample  points  each)  for  the  three  different 
house  sizes. 

Elevation  Correction 

An  effort  was  made  to  improve  the  accuracy 
of  the  weighting  factor  and  interpolation  meth- 
ods by  considering  elevation  in  the  correlations. 
It  was  hypothesized  that  most  of  the  unbiased 
variation  in  annual  heating  energy  require- 
ments not  explainable  by  latitude  and  longi- 
tude can  be  accounted  for  by  elevation.  The 
elevation  correction  for  each  of  these  methods 
was  based  on  a  second  order  multiple  regres- 
sion, using  all  24  locations  with  known  weather 
data  to  solve  for  the  7  unknown  coefficients. 
The  regression  equation  took  the  form: 

E  =  aL  +  b®  +  cH>  +  dL2  +  e4>2  +fH  +  g   (9) 


H  =  elevation,  ft. 

The  coefficient,  /,  was  the  linear  correction 
factor  accounting  for  temperature  changes  due 
to  elevation,  similar  to  the  adiabatic  lapse  rate. 
The  elevation  correction  established  by  this 
correlation  of  all  24  data  points  was  applied  to 
the  inverse  square  distance  method  and  to  each 
of  the  interpolation  methods  by  rearranging 
the  appropriate  equation  (6-9)  to  give, 


E-fH=... 


(10) 


and  solving  for  the  coefficients  as  before.  Add- 
ing back  in  the  elevation  correction  for  each  lo- 
cation gave  the  estimated  heating  energy 
requirement  corrected  for  elevation. 

DISCUSSION  OF  RESULTS 

Table  5  shows  the  standard  deviations  of  the 
heating  energy  requirement  predictions  com- 
pared to  BLAST  predictions  for  each  of  the 
methods  used.  All  standard  deviations  were 
based  on  the  sum  of  the  squared  errors  (or  per- 
cent errors)  between  interpolation  and  BLAST 
predictions  for  the  three  houses  in  all  24  test 
sites  (a  total  of  72  data  points).  Figures  3  to  8 
show  plots  of  these  predicted  results  against 
BLAST  simulation  results  for  the  significant 
methods  discussed. 

The  first  observation  was  that  the  degree- 
day  method  was  less  accurate  than  any  of  the 
other  methods  compared.  An  important  second 
observation  was  that  the  use  of  the  linear  cor- 
rection for  elevation  produced  a  distinct  im- 
provement in  most  of  the  methods  to  which  it 
was  added. 

Considering  all  the  methods  used,  it  was  ap- 
parent that  the  choice  of  method  depended  on 
the  number  of  locations  with  weather  data  for 
use  in  interpolation.  If  only  three  neighboring 
sites  with  weather  data  were  available,  then 
linear  interpolation  with  elevation  correction 
was  best.  If  six  or  seven  sites  with  weather 
data  were  available,  then  linear  interpolation 
with  elevation  correction  was  still  best  even 
though  second  order  interpolation  was  possible. 
If  many  weather  sites  (in  this  case,  24)  were 
available  then  the  second-order  interpolation 
method  with  elevation  correction  was  best. 


where, 
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Table  5.  Standard  Deviations(SD)  of  Predictive  Methods  from  BLAST  Simulation  Data 


Heating  Energy  Requirement  Prediction  Method 


SD(error) 


SD(%  error)         Rank 


Degree-Day  Method 

Inverse  Distance  Weighting  (3  sample  pts) 

Inverse  Square  Distance  Weighting  (3  sample  pts) 

Inverse  Square  Distance  Weighting  +  Elevation  Correction 

Linear  Interpolation  (3  sample  pts) 

Linear  Interpolation  +  Elevation  Correction  (3  sample  pts) 

Linear  Interpolation  (7  sample  pts) 

Linear  Interpolation  +  Elevation  Correction  (7  sample  pts) 

2nd  Order  Interpolation  (6  sample  pts) 

2nd  Order  Interpolation  +  Elevation  Correction  (6  sample  pts) 

2nd  Order  Interpolation  (7  sample  pts) 

2nd  Order  Interpolation  +  Elevation  Correction  (7  sample  pts) 

Linear  Interpolation  (24  sample  pts) 

Linear  Interpolation  +  Elevation  Correction  (24  sample  pts) 

2nd  Order  Interpolation  (24  sample  pts) 

2nd  Order  Interpolation  +  Elevation  Correction  (24  sample  pts) 
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11.3% 

16 
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7.5% 
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7.1% 
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8 
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15 
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6.6% 

6 
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7.1% 

9 
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3 
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6.7% 

7 
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8.9% 

14 
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8.1% 

12 
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8.9% 
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6.1% 
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The  standard  deviations  of  the  heating  en- 
ergy requirement  predictions,  based  on  interpo- 
lation from  BLAST  simulations  at  the  same 
sites,  ranged  from  four  to  seven  percent  for 
these  best  methods,  depending  on  the  number 
of  neighboring  sites  with  weather  data  used  in 
the  interpolations.  These  are  minor  errors  com- 
pared with  those  introduced  by  imprecise 
house  description  and  occupant  preferences 
(thermostat  settings,  natural  ventilation,  appli- 
ance usage,  etc.).  Furthermore,  these  errors 
will  be  reduced  in  real  applications  because  the 
interpolation  sites  used  in  this  study  will  be 
available  as  reference  locations  which  will 
move  the  set  of  reference  sites  closer  together. 

CONCLUSIONS 

For  design  or  research  applications  where 
building  energy  analysis  computer  programs 
are  used,  interpolation  methods  are  superior  to 
simplified  heating  energy  calculation  methods 
such  as  the  degree-day  method  in  locations 
where  detailed  weather  data  are  not  available. 
Performing  multiple  building  energy  analyses 
and  applying  the  interpolation  techniques  de- 
scribed in  this  paper  might  seem  time  consum- 
ing. However,  the  primary  task  in  building 


energy  analysis  is  still  the  production  of  input 
data  describing  the  house  to  the  computer  pro- 
gram. Running  the  program  at  multiple 
weather  sites  and  then  interpolating  between 
them  could  easily  be  automated  to  produce 
heating  energy  use  predictions  at  sites  without 
weather  data.  Standard  errors  of  four  to  seven 
percent  are  expected  for  the  better  interpola- 
tion methods,  compared  with  over  eleven  per- 
cent for  the  degree-day  method. 
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Microbial  Contamination  of  Building  HVAC  Systems: 
Causes  And  Solutions 


Ossama  A.  Abdou  and  Francis  A.  Sando 


ABSTRACT 


Modern  buildings  today  are  experiencing  in- 
door air  quality  (IAQ)  problems.  There  are 
many  causes  of  poor  air  quality,  some  of  which 
are  leading  to  a  variety  of  building-related  com- 
plaints. One  of  the  main  causes  is  indoor  micro- 
bial contamination — bioaerosols — disseminated 
throughout  the  building  by  the  air-conditioning 
system.  Bioaerosols  consist  of  viruses,  bacteria 
and  fungi,  that  thrive  in  moist  areas.  The  air- 
conditioning  system,  with  wet  cooling  coils  and 
moist  ducts,  make  ideal  host  reservoirs  for 
these  microrgamsms.  The  most  publicized  bac- 
terium associated  with  air-conditioning  sys- 
tems is  Legionella  pneumophila  which  causes 
Legionnaires'  disease  and  Pontiac  fever.  Le- 
gionella strains  have  been  found  in  water  cool- 
ing towers  associated  with  air-conditioning 
plants  and  in  water  spray  systems.  Sapro- 
phytic organisms  such  as  Aspergillus  flavus 
thrive  in  some  nonliving  plant  substrates  and 
can  be  pulled  into  a  building  through  the  fresh 
air  intake  grille  of  the  air-conditioning  system. 
They  are  known  to  have  a  carcinogenic  effect. 
The  extent  of  microbial  contamination  in  build- 
ing HVAC  systems  is  described  and  remedial 
actions  are  proposed. 
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Today's  air-conditioning  system  can  no 
longer  supply  only  heated  or  cooled  air  to  pro- 
vide occupant  comfort,  but  also  must  provide 
systems  to  preserve  the  quality  of  the  air  we 
breathe.  The  system  must  have  components  to 
kill  bioaerosols,  as  well  as  eliminate  chemical 
and  particulate  contamination.  A  generic  air- 
conditioning  system,  with  such  components 
shown,  is  proposed  and  diagramed  in  this  pa- 
per. 

INTRODUCTION 

Indoor  air  pollution  is  an  emerging  problem 
which  has  been  linked  to  energy  conservation, 
the  use  of  building  materials  and  the  design 
and  operation  of  heating,  ventilation  and  air- 
conditioning  (HVAC)  systems.  Since  a  large 
amount  of  time — up  to  90  percent — is  spent  in- 
doors (EPA  1989;  Chen  et  al.  1992),  indoor  air 
quality  (IAQ)  is  an  important  factor  to  our 
health  and  comfort.  IAQ  problems  and  com- 
plaints are  being  increasingly  reported  (ACGIH 
1984;  ORNL  1985;  Hodgson  and  Hess  1992) 
and  have  led  to  the  classification  of  the  follow- 
ing building  sicknesses:  sick  building  syndrome 
(SBS),  building-related  illness  (BRI),  and  multi- 
ple-chemical sensitivity  (MCS).  Increased 
awareness  of  the  potential  health  risks  associ- 
ated with  indoor  air  contaminants  arguably 
linked  with  significant  lost  annual  worker  pro- 
ductivity along  with  related  worker  claims  and 
litigations,  has  stimulated  the  need  to  improve 
our  understanding  of  how  contaminants  are 
transported  in  buildings.  IAQ  is  a  result  of  a 
complex  relationship  between  the  contamina- 
tion sources  in  a  building,  the  ventilation  rate, 
sensitivity  of  various  individuals  to  particular 
contaminants,  and  the  dilution  of  the  indoor  air 
contaminant  concentrations  with  outdoor  air  or 
treated  recirculated  air.  This  complex  relation- 
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ship  is  further  aggravated  by  outdoor  sources 
used  for  dilution  air  and  pollution  sinks  in  a 
building  that  may  amplify  or  act  as  a  host  for 
contaminants.  One  important  implication  of 
this  is  that  building  HVAC  systems  which  cre- 
ate and  control  interior  environments  can  pro- 
mote or  prevent  and  mitigate  indoor  air 
contamination  through  their  design,  operation 
and  maintenance. 

CAUSES  OF  BUILDING  SICKNESS 

Firstly,  the  term  SBS  is  a  misnomer,  be- 
cause the  buildings  are  not  sick,  although  occu- 
pants of  such  buildings  exhibit  symptoms 
associated  with  being  ill.  The  World  Health  Or- 
ganization identifies  a  range  of  symptoms  for 
SBS  which  cause  distress  to  some  building  occu- 
pants, but  cannot  be  clinically  diagnosed  and, 
therefore,  cannot  be  medically  treated.  Typi- 
cally, in  such  a  case,  a  building's  occupants  ex- 
hibit discomfort  or  dissatisfaction  with  the 
environment  other  than  feeling  too  warm  or  too 
cold.  Occupant  discomfort  is  usually  evidenced 
by  symptoms  such  as  headache,  fatigue,  sinus 
congestion,  sore  throats,  and  eye  irritation.  The 
occupants  experience  relief  almost  immediately 
after  they  exit  the  building  and  the  symptoms 
disappear  when  the  sufferers  leave  the  building 
for  the  weekend.  Usually,  the  symptoms  recur 
upon  re-entry. 

This  definition  of  SBS  implies  that  sensory 
irritation  must  be  a  dominant  complaint  while 
systemic  symptoms  should  be  infrequent.  It  fur- 
ther implies  that  the  prevalence  of  the  symp- 
toms should  be  so  great  that  the  reactions  are 
clearly  related  to  an  unfavorable  indoor  envi- 
ronment and  not  caused  by  any  medical  condi- 
tions among  the  occupants.  Unless  the 
dissatisfaction  level  is  around  20  percent  of  the 
occupants,  the  problems  are  not  likely  to  be  con- 
sidered as  SBS,  and  if  the  percentage  dissatis- 
fied rises  above  20  percent  in  one  part  of  the 
building,  where  the  environment  is  said  to  be 
common  throughout,  the  problem  is  probably 
much  more  mundane. 

BRI,  on  the  other  hand,  refers  to  an  illness 
caused  as  a  result  of  exposure  to  an  indoor  air 
agent  that  can  be  identified.  As  it  is  commonly 
used,  the  term  BRI  refers  to  an  illness  out- 
break among  building  occupants  (e.g.,  individu- 
als entering  a  building  for  extended  periods, 


and  those  residing  and  working  in  a  particular 
building),  with  no  secondary  spread  of  the  ill- 
ness to  others  outside  the  building  with  whom 
affected  individuals  come  into  contact.  Symp- 
toms include  infection,  hypersensitivity  pneu- 
monitis (HP),  etc.  MCS  refers  to  high  sensitiv- 
ity to  chemically-based  pollutants  in  the  envi- 
ronment. 

Building  sicknesses  (especially  SBS)  have 
been  initially  ascribed  to  inadequate  ventila- 
tion— either  because  of  increased  envelope 
tightness  or  decreased  introduction  of  outdoor 
air  into  a  building  space.  A  number  of  other  fac- 
tors have  been  identified  as  causing  sickness  of 
occupants.  These  include  exposure  to 
bioaerosols,  volatile  organic  compounds 
(VOCs),  unpleasant  odors,  smoke,  chemicals 
and  particulates  found  in  the  air  of  both  natu- 
rally and  mechanically  ventilated  buildings. 
One  class  of  indoor  air  pollutants  are  biological 
contaminants  that  thrive  in  moist  areas,  includ- 
ing air-conditioning  systems.  Major  causative 
agents  contributing  to  a  large  number  of  "sick 
buildings"  have  been  reported  as  being  primar- 
ily fungal  and  bacterial  contamination  (Cun- 
ningham 1989).  Fungi  generate  not  only 
spores,  but  also  release  toxic  compounds  and 
VOCs  in  the  air.  Since  these  organisms  are  be- 
ing implicated  in  an  increasing  number  of  prob- 
lem buildings,  not  only  a  better  understanding 
of  their  nature,  habitat  and  control  is  war- 
ranted, but  also  remedial  actions  are  to  be 
sought  to  mitigate  and  eventually  eradicate  the 
problem  of  microbial  contamination  in  building 
HVAC  systems. 

It  is  important  to  recognize  that  microbial 
contamination  can  occur  in  a  building  through 
sources  other  than  the  HVAC  system.  For  ex- 
ample, flooding,  roof  leaks  and  some  janitorial 
practices  (e.g.,  carpet  shampooing)  have  been 
referred  to  as  contributing  to  such  microbiologi- 
cal growth  in  building  spaces.  The  existence  of 
microbial  indoor  air  pollutants  is  almost  al- 
ways evidence  that  the  building's  maintenance 
or  repair  procedures  are  lacking  in  some  re- 
spect. 

MICROBIAL  CONTAMINATION 

Microbial  contamination  is  a  catch-all  term 
for  any  bioaerosol  consisting  of  airborne  parti- 
cles, viruses,  bacteria,  fungi,  oversized  mole- 
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cules  or  volatiles  which  are  living  or  released 
into  the  air  by  some  living  organism. 
Bioaerosols  include  vegetative  microbial  cells, 
their  reproductive  units  and  metabolites  that 
are  small  enough  or  sufficiently  volatile  to 
achieve  aerial  dispersion.  Bioaerosols  range 
widely  in  size,  from  less  than  0.1  micron  to 
greater  than  100  microns.  Three  conditions  are 
usually  necessary  to  release  them  into  the  air, 
a  reservoir,  amplification  or  a  means  to  in- 
crease their  concentration,  and  dissemination. 
Domestic  animals,  birds,  and  humans  emit  or- 
ganisms or  effluents  which  can  easily  act  as  res- 
ervoirs and  amplifiers.  The  other  conditions, 
amplification  and  dissemination,  occur  from 
the  living  host.  Influenza,  a  respiratory  infec- 
tion, is  typically  spread  in  this  manner. 

Levels  of  bacteria  and  fungi  in  air  are  meas- 
ured in  terms  of  colony  forming  units  (CFUVm3 
or  in  concentrations  of  magnitude/mL.  A  colony 
typically  consists  of  106  to  107  cells.  Magni- 
tudes of  100  CFU/m3  or  less  for  some  bacteria 
are  considered  low.  Indoor  concentrations  of  as 
high  as  100,000  CFU/m3  are  not  uncommon  in 
large  buildings  (Morey  1988a).  The  most  com- 
mon outdoor  composition  of  airborne  fungi  are 
as  follows:  Cladosporium  sp.,  Penicillium  sp., 
Alternaria  sp.  and  Aspergillus  sp.  Other  pre- 
dominant microorganisms  include  Rhodotorula 
sp.,  Legionella  pneumophila,  Aureobasidium 
sp.,  Acanthamoeba  polyphaga  and  Thermoacti- 
nomyces  (Hodgson  et  al.  1985).  The  transport 
mechanism  for  viruses  is  by  means  of  droplets 
(e.g.,  generated  by  sneezing  or  coughing)  or 
clinging  to  dust  particles. 

The  most  publicized  bacterium  associated 
with  air-conditioning  systems  is  Legionella 
pneumophila,  which  causes  Legionnaires'  dis- 
ease and  Pontiac  fever.  More  than  25  species  of 
this  bacterium  have  now  been  recognized  with 
many  more  serotypes.  The  bacterium  occurs  in 
nature  in  moist  soil  or  tepid  water — up  to  50 
percent  of  the  bodies  of  water  in  the  US  are  be- 
lieved to  contain  Legionella  strains  (Hodgson 
and  Hess  1992) — and  has  been  found  in  water 
cooling  towers  associated  with  air-conditioning 
plants  and  in  water  spray  systems  which  form 
part  of  some  air-conditioning  equipment.  Le- 
gionella must  have  a  water  substrate  to  live, 
and  it  can  be  dangerous  if  water  droplets  are 


picked  up  by  the  HVAC  system  and  carried 
into  an  occupied  space.  The  level  of  risk  to  occu- 
pants is  still  considered  to  be  small,  although  a 
number  of  fatalities  have  been  reported.  Cer- 
tain conditions  for  the  growth  of  Legionella 
pneumophila  are  necessary,  and  although  they 
can  be  found  in  bodies  of  water  ranging  in  tem- 
perature up  to  140°F  (60°C),  significant  multi- 
plication is  generally  restricted  to 
temperatures  of  68°F  (20°C)  to  115°F  (46°C). 
Some  investigators  have  proposed  a  practical 
working  level  of  10,000  Legionella  pneumo- 
phila per  milliliter,  above  which  serious  consid- 
eration should  be  given  to  cleaning  and 
chemical  disinfection  (Fliermans  and  Nygren 
1987). 

Microbial  sampling  undertaken  in  a  number 
of  studies  has  shown  that  the  concentration  of 
fungi  present  in  indoor  air  can  be  higher  than 
that  found  outdoors  by  an  order  of  magnitude 
or  more.  This,  in  turn,  suggests  that  internal 
sources  of  fungi  exist  and  are  commonly  attrib- 
uted to  incursion  of  moisture  into  various  com- 
ponents within  the  space  (e.g.,  carpets,  ceiling 
tiles,  upholstered  furniture,  wall  partitions  and 
office  materials)  or  to  the  HVAC  system  itself 
as  described  below.  A  few  large-scale  studies 
have  been  published  and  have  demonstrated 
that  the  complaints  and  symptoms  are  more 
frequent  in  buildings  with  mechanical  ventila- 
tion (Mendell  and  Smith  1990).  On  the  other 
hand,  mechanical  ventilation  is  usually  associ- 
ated with  a  lowering  of  the  indoor  concentra- 
tion of  airborne  fungi  compared  with  levels 
encountered  outdoors  or  in  naturally  ventilated 
spaces  (Hirsch  et  al  1978;  Rose  and  Hirsch 
1979;  Solomon  et  al.  1980).  Table  1  lists  a  few  of 
the  fungi  which  can  affect  building  occupants 
by  collecting  near  fresh  air  intakes,  at  filters, 
near  cooling  coils  or  within  the  building  occu- 
pied space. 

Because  of  the  many  bioaerosols  that  can 
contaminate  the  air  we  breathe  and  make  us 
sick,  the  air  handling  units  (AHU)  and  duct- 
work of  a  building  air-conditioning  system  are 
now  coming  under  scrutiny.  Is  the  air-condition- 
ing system  a  reservoir  or  an  amplifier?  It  cer- 
tainly could  be  a  disseminator  since  it  is 
delivering  air  to  (and  collecting  air  from)  occu- 
pied spaces. 
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Table  1 .  Sources  and  Health  Efecls  of  Fungi  Commonly  Found  in  Buildings 


Organism 


Causal  Part 


Main  Sources 


Health  Effects 


Aspergillus 


Spores 


Damp  organic  material; 
Self-heated  compost 


Invasive  aspergillosis; 
Aspergilloma;  Allergic 
bronchopulmonary 
aspergillosis  (ABPA) 


Aspergillus  Penicillium 
Sporobolomyces 

Spores 

Damp  organic  material; 
Standing  water 

Hypersensitivity 
pneumonitis  (HP) 

All  airborne  fungi 
spores 

Spores 

Damo  organic  material; 
Outdoor  air 

Allergic  asthma; 
Rhinitis 

Aspergillus  flavus 

Toxin 

Agricultural  products 

Cancer 

Cryptococcus 

Spores 

Bird  droppings 

Cryptococcosis 

Histoplasma 

Spores 

Bird  droppings 

Histoplasmosis 

Penicillium  species 

VOC 

Damp  organic  material 

Irritants 

Stachybotrys  atra 

Toxin 

Damp  cellulosic  material 

Acute  Toxicosis 

Ulocladium 

VOC 

Damp  organic  material 

Headache 

MICROBIAL  PROLIFERATION  THROUGH 
HVAC  SYSTEMS 

Microbial  contaminants  usually  are  brought 
from  their  sources  to  the  occupied  space  by 
means  of  an  airstream.  Inside  a  building,  con- 
taminants are  distributed  from  their  sources  by 
means  of  natural  convection  or  the  ventilation 
system  itself.  Many  studies  show  that  ventila- 
tion systems  play  a  major  role  in  IAQ  prob- 
lems. Over  50  percent  of  all  IAQ  problems 
investigated  over  an  18  year  period  by  the  Na- 
tional Institute  for  Occupational  Safety  and 
Health  (NIOSH)  have  been  attributed  to  venti- 
lation deficiencies  (Moseley  1990).  Regulation 
of  fresh  air  quantities  and  air  movement  is  sug- 
gested as  a  measure  in  alleviating  most  reoccur- 
rences of  SBS.  It  also  appears  that  air  flow 
patterns  in  a  building  may  have  a  significant  in- 
fluence on  the  spreading  of  airborne  bacteria 
and  viruses.  This  emphasizes  the  importance  of 
ventilation  effectiveness  (whether  mechanical 
or  natural)  in  buildings.  On  the  other  hand, 
since  the  primary  requirement  for  the  growth 
and  amplification  of  microorganisms  in  the  in- 
door environment  is  the  presence  of  moisture, 
contamination  can  occur  in  a  building  regard- 
less of  the  amount  of  outdoor  air  being  supplied 
to  the  space.  The  exact  number  of  problem 
buildings  with  microbial  contamination  is  not 
known;  however,  it  is  estimated  that  these  can 


be  as  high  as  33  percent  (Wallingford  1986) 
(HBI  1990). 

Microbial  contaminants  in  the  indoor  envi- 
ronment may  cause  illnesses  of  two  general 
types,  namely,  those  that  are  infective  and 
those  that  are  allergic  in  nature.  Many  of  these 
illnesses  occur  as  a  response  to  antigens 
aerosolized  from  HVAC  system  components 
such  as  humidifiers  or  air  washers  which  may 
be  contaminated  (Edwards  1980;  Scully  et  al. 
1979).  Other  sources  of  reported  illnesses  are 
from  other  building  components  that  have  been 
damaged  by  recurrent  floods  or  other  sources  of 
moisture  (Morey  1988a). 

The  presence  of  microbes  in  the  indoor  envi- 
ronment in  sufficient  amounts  or  types  to  cause 
illnesses  depends  on  a  number  of  factors: 

•  type  of  mechanical  components 
that  constitute  the  HVAC  system; 

•  maintenance  plan  for  mechanical 
systems;  and 

•  the  presence  of  water. 

Water  is  the  primary  factor  permitting  mi- 
crobial contamination.  Since  absence  of  water 
or  moisture  is  virtually  impossible  in  any 
HVAC  system,  the  AHUs  of  such  systems  offer 
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ideal  environments  for  microbial  proliferation. 
HVAC  systems  must  function  to  provide  an  en- 
vironment where  temperature,  air  distribution, 
humidity  and  IAQ  are  maintained  within  the 
described  levels.  For  the  most  part,  HVAC  sys- 
tems are  designed  to  mix  a  given  amount  of  out- 
door air  with  another  amount  of  recirculated 
air,  condition  this  mixture,  and  distribute  it  to 
the  occupied  space.  By  mixing  outdoor  air  with 
recirculated  air,  dilution  of  pathogenic  bacteria 
is  possible  as  well  as  reduction  of  chemical  vi- 
tiation. 

Outdoor  air  entering  the  building  may  be 
contaminated  by  nearby  sources  of  bioaerosols 
such  as  those  from  cooling  towers,  evaporative 
condensers,  toilet  exhausts  and  sanitary  vents. 
Cooling  system  components  are  especially 
known  to  harbor  microorganisms  which  can 
cause  a  variety  of  problems,  including  loss  of 
heat  transfer  efficiency  and  deterioration  of  sys- 
tem components.  They  have  further  been  impli- 
cated in  the  growth  and  dissemination  of 
pathogenic  microorganisms.  Thus,  a  HVAC  sys- 
tem may  serve  to  transport  microorganisms 
from  the  locus  of  contamination  to  the  vicinity 
of  sensitive  occupants.  The  air  exiting  the  mix- 
ing plenum  of  the  AHUs  passes  through  filters. 
These  filters,  depending  on  their  efficiency,  are 
poor  collectors  of  small  spores  that  might  be  en- 
tering the  AHU  from  an  outdoor  source  such  as 


dead  vegetation.  The  filtered  mixture  of  out- 
door and  return  air  passes  through  the  coil  sec- 
tion where  it  is  heated  or  cooled  to  meet  the 
specified  requirements.  Conditioned  air  is  then 
supplied  to  the  building  spaces.  This  is  shown 
conceptually  in  Figure  1 .  A  schematic  diagram 
of  a  blow-through  air  handler  in  which  the  fan 
is  located  upstream  of  the  heating  and  cooling 
coils  is  given  in  Figure  2. 

During  the  summer  season,  moisture  con- 
denses from  the  air  stream  as  it  passes  through 
the  cooling  coils  [which  usually  maintain  a  42- 
55°F  (6-13°C)  chilled  water  system  tempera- 
ture or  a  30  to  45°F  (-1  to  7°C)  surface 
temperature  of  a  DX  system].  When  the  coil 
surface  temperature  is  below  the  dewpoint  of 
the  air  stream,  condensed  water  collects  in  a 
drain  pan  and  exits  the  AHU  through  drain 
lines.  In  many  cases,  water  stagnates  in  the 
drain  pans;  hence,  the  pan  becomes  a  potential 
reservoir  and  amplification  site  for  microbes. 
Furthermore,  air  washer  type  humidifiers  lo- 
cated downstream  from  filter  banks  and  near 
heating  and  cooling  coils  operate  by  aerosoliz- 
ing water  into  the  air  stream.  Water  that  is  not 
taken  up  by  the  air  stream  falls  into  a  sump 
and  is  circulated  to  nozzles  for  re-aerosoliza- 
tion.  This  type  of  humidification  system  is  par- 
ticularly subject  to  extensive  microbial 
contamination  as  well  as  being  an  ineffective 
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Figure  1 .  Concept  of  a  typical  HVAC  system 
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Figure  2.  Schematic  of  a  blow-through  air  handler 
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Figure  3.  Conventional  spray-type  air  washer 


way  to  humidify  air.  In  Figure  3  a  typical  air 
washer  is  illustrated. 

Conditioned  air  leaving  the  AHU  fan  plenum 
is  transported  to  the  space  by  a  system  of 
ducts.  Conditioned  air  from  the  AHU  enters 


the  space  through  grilles  or  diffusers  and  mixes 
with  room  air.  Mixed  air  is  then  pulled  out  of 
the  space  either  through  a  return  duct  or 
through  a  common  return  plenum  and  is  then 
partially  removed  by  the  return  fan  and  is  par- 
tially exhausted.  Because  of  difficulty  in  access- 


Volume  3,  Number  1,  1994 


2b' 


Fin  tube  coil 


Metal  enclosure 


Condensate 
drain  pan 


Air  Blower 
and  motor 


Return  air  damper 

with  filter 
Drain  line 

Figure  4.  Typical  fan  coil  unit  (below  window  type) 


ing  the  return  plenum  and  ducts  for  proper 
maintenance,  microbial  contamination  that 
may  occur  is  removed  only  with  great  difficulty. 
Ducts  that  have  been  installed  and  operating 
for  ten  to  twenty  years  without  cleaning  are 
probably  loaded  with  dust  and  dust  mites.  Ceil- 
ing plenums  are  particularly  susceptible  to 
dust  and  fibers  which  are  then  pulled  into  the 
return  duct  connecting  to  the  AHU. 

The  heating  and  cooling  load  of  a  large  build- 
ing may  be  handled  by  fan  coil  units  installed 
at  the  periphery  of  the  building.  Each  unit  con- 
sists of  low-efficiency  filters,  a  fan,  a  heating 
and  cooling  coil  and  a  drain  pan.  Water  may  ac- 
cumulate in  the  pan  and  create  a  microbial 
slime.  Figure  4  shows  a  typical  fan  coil  unit  de- 
picting the  location  of  the  drain  pan.  Porous  in- 
sulation liners  applied  to  the  inside  of  the  fan 
coil  units  for  sound  attenuation  provide  an  ad- 
ditional substrate  for  microbial  growth.  Dirt 
and  debris  accumulate  within  the  insulation 
and  microorganisms  may  flourish  on  this  sub- 
strate when  it  becomes  moist  during  the  air- 
conditioning  mode  of  operation.  If  mold  or 
slime  is  visible,  then  fungi  are  present. 

INDOOR  MICROBIAL  EXPOSURE  LEVEL 

No  scientifically  based  criteria  exist  to  show 
whether  a  measured  level  of  fungi  or  bacteria 
is  a  risk  factor  with  regard  to  any  allergic  dis- 
eases. Therefore,  a  carefully  thought-out  and 
planned  approach  must  be  used  to  determine 
these  risks.  Any  quantitative  criteria  must 


take  into  account  the  qualitative  nature  of  the 
diverse  viable  and  nonviable  etiologic  agents 
thought  to  be  responsible  for  these  illnesses. 
Several  suggestions  have  been  made  concern- 
ing normal  or  hygienically  acceptable  levels  of 
airborne  viable  particulates  found  in  large 
buildings  and,  in  particular,  large  office  envi- 
ronments. A  level  of  approximately  1,775  bacte- 
ria-containing particles/m3  was  described  as 
the  threshold  for  clerical  environments  in  need 
of  investigation  and  hygienic  improvement 
(Bourdillon  et  al.  1948).  Others  suggested  that 
this  threshold  be  1000  CFU/m3  (Morey  et  al. 
1984).  Levels  of  about  700  bacteria/m3  were  con- 
sidered more  normal,  and  it  has  been  stated 
that  levels  of  viable  microorganisms  in  the  in- 
door environment  seldom  exceed  1,700/m3 
(Wright  et  al.  1969).  Other  studies  indicate 
that  levels  of  viable  particulates  in  office  spaces 
with  HVAC  systems  normally  range  from  100 
to  300  CFU/m3  (Yoshizawa  and  Sugawara 
1985).  The  ACGIH  (1989)  has  developed  an  ap- 
proach that  bases  the  allowable  indoor  contami- 
nant level  on  existing  outdoor  levels. 

In  reporting  on  a  study  on  fungi  in  Canadian 
houses,  Miller  et  al.  (1988)  proposed  the  follow- 
ing criteria  of  acceptability  for  fungal  pollution 
in  indoor  air:  More  than  50  CFU/m3  should  be 
reason  for  concern  if  there  is  only  one  species 
present,  while  less  than  150  CFU/m3  should  be 
considered  acceptable  if  there  is  a  mixture  of 
species,  and  less  than  300  CFU/m3  should  be  ac- 
ceptable if  the  species  present  are  primarily 
Cladosporium  or  other  common  phylloplane 
fungi.  More  recently,  a  study  recommended 


that  the  airborne  microbe  levels  be  between 
250  and  750  CFU/m3  (HBI  1991)  to  avoid  any 
health  hazards.  Another  study  indicates  that  a 
factor  denning  indoor  contamination  includes 
indoor  bioaerosol  levels  being  consistently  more 
than  double  levels  outdoors  and  exceeding  1000 
spores/m3  of  air  (Burge  1990).  Nonetheless, 
many  studies  point  to  the  inappropriateness  of 
using  any  given  level  or  concentration  of  viable 
particulates  as  the  sole  indicator  of  an  illness 
threshold  (Morey  1988a). 

A  Norwegian  Health  Directorate  Report 
(1990)  indicated  a  very  general  or  rather  vague 
threshold,  characterizing  exposure  limits  for 
microorganisms.  It  simply  stated  that  a  mouldy 
smell  should  not  occur.  There  is  no  question 
that  using  the  sense  of  smell  as  a  tool  to  iden- 
tify maximum  contaminant  limits  is  a  practice 
fraught  with  inaccuracies  as  well  as  uncertain- 
ties. 

It  is  indeed  interesting  to  note  that  despite 
much  research,  the  direct  relationship  between 
many  pollutants  and  health  problems  cannot 
be  proven.  However,  there  is  sufficient  evi- 
dence to  warrant  some  action  in  controlling  the 
levels  of  microbial  contaminants.  For  example, 
it  has  been  reported  that  the  levels  of  contami- 
nation posed  a  danger  in  27  percent  of  "sick 
buildings"  investigated  as  a  result  of  excessive 
fungal  contamination  and  in  25  percent  as  a  re- 
sult of  excessive  bacterial  contamination  (Cun- 
ningham 1989).  There  is  some  correlation 
between  concentration  levels  of  microbial  pol- 
lutants and  mortality  and  illness  rates.  Never- 
theless, no  single  pollutant  has  been  singled 
out  as  being  responsible.  Further  research  is 
proceeding  to  quantify  the  risks  and  identify 
those  parameters  which  influence  risk.  Since 
most  microbial  contaminants  affect  the  human 
body  by  means  of  the  respiratory  system  and 
most  fungi  are  facultative  parasites,  there  is  a 
close  link  between  the  HVAC  system,  air  qual- 
ity and  complaints  about  SBS.  It  is  suggested 
that  there  is  a  synergistic  association  among  a 
number  of  factors  including  multiple  exposure 
levels  to  pollutants  and  an  individual's  general 
state  of  health. 

MICROBIAL  CONTROL  IN  HVAC  SYSTEMS 

Microbes  are  ubiquitous  in  nature  and  are 
normally  found  on  environmental  surfaces. 


They  enter  a  building  in  many  ways,  and  in- 
deed may  be  built  into  it,  because  construction 
dust  and  debris  may  contaminate  HVAC  sys- 
tem components  at  the  time  of  their  installa- 
tion. Since  microbes  and  their  spores  are  small 
and  light  enough  to  be  carried  with  air  cur- 
rents, the  outside  air  taken  into  a  building  al- 
ways contains  large  numbers  of  them.  The 
outdoor  air  intake  of  an  HVAC  system  should, 
therefore,  be  protected  with  a  filter  or  a  grille 
of  fine  mesh  to  prevent  organic  materials  from 
being  drawn  in,  these  materials  being  particu- 
larly potent  in  microbial  contamination. 

The  ventilation  system  in  the  building  pro- 
vides the  first  line  of  defense  in  arresting  mi- 
crobes. The  return  air,  drawn  back  to  the  air 
handling  unit,  may  be  heavily  laden  with  air- 
borne particulates  and  microbes.  A  properly  de- 
signed and  maintained  system  will  have  the 
return  air  diluted  with  outdoor  air  and  both 
will  be  efficiently  filtered  to  reduce  the  microbe 
content.  In  many  cases,  however,  the  filtration 
systems  are  poorly  designed  and  installed  and, 
above  all,  badly  maintained.  This  allows  micro- 
organisms (especially  fungi)  to  amplify  and  con- 
taminate the  air  supply  system.  Standing 
water  around  cooling  coils,  condensate  trays, 
humidifiers  and  condensation  water  in  ducts, 
fan  chambers  and  on  insulation  surfaces,  then 
become  contaminated  and  allow  further  prolif- 
eration of  growth.  The  ventilating  system  itself 
then  provides  a  means  of  transporting  more 
spores  into  the  occupied  areas.  Based  on  the 
above,  prevention  from  infestation  is  of  para- 
mount importance.  Many  "spray-on"  antimicro- 
bials are  currently  in  use;  few  of  them, 
however,  give  durable  protection  since  spray-on 
compounds  wear  off  with  use,  or  they  become 
ineffective  when  cleaned  or  exposed  to  mois- 
ture. Moreover,  sprays  themselves  use  aerosol 
propellants  or  contain  components  which  in- 
crease the  chemical  burden  of  the  indoor  air 
with  potentially  dangerous  VOCs. 

Several  researchers  have  investigated  the  ef- 
ficacy of  biocides  in  controlling  microorganisms 
in  HVAC  system  components,  especially  in  cool- 
ing tower  environments.  For  example,  a  study 
demonstrated  that  a  number  of  biocides  kill  a 
variety  of  microbes  habitating  building  cooling 
systems  including  Legionella  when  tested  in  a 
laboratory  environment  (Soracco  et  al.  1983). 
Another  study  showed,  however,  that  some  of 
these  same  biocides  are  not  as  effective  in  real 
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world  situations  (Soracco  and  Pope  1983),  espe- 
cially when  concentrations  of  these  chemicals 
are  reduced  in  the  water.  Nonetheless,  it  is  ar- 
gued that  continuous  treatment  with  biocides 
that  are  effective  at  low  levels  may  be  a  viable 
treatment  option  (Pope  and  Dziewulski  1992). 
However,  the  economics  of  a  continuous  bioci- 
dal  treatment  have  not  yet  been  established. 

A  study  demonstrated  that  ozone,  if  properly 
applied,  could  affect  control  of  microbial  popula- 
tions, including  Legionella,  in  a  variety  of  cool- 
ing systems  (Pope  et  al.  1984).  Of  course,  the 
degree  to  which  microbial  communities  in  the 
HVAC  system  can,  or  should  be,  controlled  is 
opinionated.  Some  investigators  are  of  the  opin- 
ion that  maintaining  a  system  in  a  totally  ster- 
ile condition  is  impractical  if  not  impossible,  let 
alone  economical.  Instead,  the  target  should 
rather  be  to  keep  the  microorganisms  under 
control  and  to  limit  their  effects  on  operation  of 
the  system.  Experience  shows  that  this  is  not 
an  easy  task  given  the  high  growth  rate  of  mi- 
crobes and  the  non-existence  of  any  standard 
threshold  limit  values  for  microbial  exposure 
levels. 

REMEDIAL  ACTIONS 

Environmental  factors  that  contribute  sig- 
nificantly to  indoor  bioaerosol  problems  include 
the  microbial  content  of  the  outdoor  air  (which 
is  never  completely  free  of  microbial  elements), 
ventilation  mode,  indoor  occupant  density  and 
ventilation  rate,  and  the  presence  of  excessive 
heat  and  moisture  in  the  indoor  environment. 
A  number  of  remedial  actions  are  available 
that  can  effectively  reduce  or  control  indoor  mi- 
crobial contamination  through  the  HVAC  sys- 
tem. 

Three  methods  of  air  quality  control  can  be 
identified:  source  control,  removal  control  and 
dilution  control  or  clean  up.  The  former  can  be 
implemented  through  isolation,  containment 
and  local  exhaust  while  dilution  control  is  a 
means  of  reducing  the  concentrations  of  con- 
taminants in  the  environment  by  exchanging 
the  room  air  with  non-contaminated  air  or 
with  air  of  less  concentration  of  the  contami- 
nant. Removal  control  is  implemented  via  air 
cleaning  devices  placed  in  the  circulating  air- 
stream  of  an  air-conditioned  or  mechanically 
ventilated  space  and  it  is  usually  the  only 


means  available  when  it  is  necessary  to  control 
the  environment  at  concentrations  below  those 
of  the  outdoor  air.  Among  the  three  control 
methods,  source  control  is  the  only  means  to 
prevent  exposure  to  a  contaminant  altogether 
while  in  the  other  control  methods  exposure  to 
the  contaminant  is  occurring  or  assumed.  In 
other  words,  only  source  control  can  eliminate 
exposure  of  the  contaminant;  the  other  two 
methods  inherently  require  mixture  of  the  con- 
taminant with  the  indoor  air  before  any  control 
methods  are  applied.  Removal  can  also  mean 
inactivating  the  air-conditioning  system  and 
physically  washing  the  system  with  a  biocide. 

Remedial  measures  apply  primarily  to  exist- 
ing systems,  and  such  systems  can  take  two 
forms: 

•  Preserve  the  inside  occupied  space 
free  from  contamination;  and 

•  Clean  up  existing  contamination. 

A  third  measure,  applying  to  new  construc- 
tion, requires  design  professionals  to  design 
buildings  and  mechanical  systems  to  prevent 
indoor  air  quality  deterioration.  Architects  and 
mechanical  engineers  must  work  together  at 
the  conceptual  stage  to  design  and  install 
HVAC  equipment  and  ductwork  in  such  a  way 
that  bacteria  and  fungi  cannot  get  a  foothold 
for  growth  within  the  mechanical  system  or  the 
inside  building  substrate. 

The  ability  of  a  ventilation  system  to  control 
the  concentrations  of  contaminants  within  ac- 
ceptable levels  at  the  exposure  site  of  the  occu- 
pants is  highly  dependent  on  the  air 
distribution  patterns  both  within  and  between 
functional  spaces.  In  some  spaces,  natural  ven- 
tilation may  be  the  primary  method  of  dilution 
control.  In  most  cases,  the  primary  method  of 
control  is  mechanical  ventilation.  In  either 
case,  the  effectiveness  of  the  system  for  air 
quality  control  is  dependent  upon  the  system 
characteristics,  the  room  air  exchange  rate  and 
the  uniformity  of  the  airflow  patterns  within 
the  room.  If  the  room  air  distribution  is  not  suf- 
ficient to  dilute  or  remove  contaminants  from 
the  location  of  most  likely  exposure,  the  effec- 
tiveness of  the  system  will  be  impaired.  Thus, 
the  air  distribution  patterns  within  the  room 
may  be  as  important  to  the  effectiveness  of  the 
ventilation  system  as  the  room  air  exchange 
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rate.  Systems  which  characteristically  provide 
little  air  movement  (some  VAV  systems)  can  ex- 
acerbate contamination  by  bioaerosols. 

Certain  heating  systems  may,  by  virtue  of 
their  operation,  be  more  inviting  to  microbial 
contamination  than  others.  For  example,  forced- 
air  heating  systems  may  allow  accumulation  of 
spores  in  ductwork  and  produce  bioaerosols  by 
subsequent  reflotation  of  such  particles,  as  well 
as  particles  contained  in  settled  dust  within  the 
living  space.  In  contrast,  radiant  heat  imposes 
neither  such  inaccessible  ductwork  nor  the 
brisk  air  movement  (i.e.,  ventilation)  presented 
by  forced-air  systems.  Heating  systems  based 
on  combustion  of  fossil  fuels  exhaust  poten- 
tially toxic  combustion  products  to  the  outside 
and  may  cause  substantial  outside  air  to  infil- 
trate into  a  building  through  cracks. 

Preservation  of  uncontaminated  space 

How  can  the  inside  space  be  kept  free  of  con- 
tamination? The  answer  is  by  applying  most  of 
the  design  parameters  engineers  are  familiar 
with  today,  plus  applying  a  few  "new"  ones 
which  are  mostly  common  sense.  Design  pa- 
rameters engineers  already  know  are  outlined 
below: 

•  Use  adequate  filtration.  Filters  50 
to  70  percent  efficient  (based  on 
ASHRAE  standards)  will  remove 
most  microorganisms  from  the  air 
stream  since  they  are  usually  up  to 
5  microns  in  diameter.  The  use  of 
upstream  inexpensive  roughing  pre- 
filters  extend  the  life  of  the  filters. 

•  Keep  the  building  under  positive 
pressure  to  reduce  contamination 
entering  through  building  envelope 
cracks  and  crevices. 

•  Use  outdoor  air  to  dilute  vitiation 
caused  by  air  recirculation.  Energy 
costs  for  heating  and  cooling  out- 
door air  can  be  significantly  lowered 
by  using  an  air-to-air  heat  ex- 
changer or  run-around  coils.  The 
use  of  "cooling  wheels"  is  question- 
able, since  microbial  contamination 
can  easily  be  transferred  along  with 
heat. 


Locate  fresh  air  intakes  remote 
from  cooling  towers,  evaporative 
condensers,  toilet  room  exhaust 
ducts,  sanitary  plumbing  vents,  in- 
dustrial exhaust  fans,  rotating  or 
gravity  exhaust  vents,  process 
equipment  vent  ducts,  pressure  re- 
lief valve  stack  heads,  pathology 
room  exhausts,  autoclave  exhausts, 
or  any  other  hospital  exhaust  ducts. 
Moreover,  the  fresh  air  intake 
should  be  located  upwind  from  all 
exhaust  ducts. 

Humidifiers  should  be  steam,  lo- 
cated downstream  from  a  pressure 
reducing  valve  to  take  advantage  of 
superheat,  and  generated  from  a 
special  small  boiler  or  a  heating 
plant  boiler  free  from  water  treat- 
ment chemicals.  The  use  of  water 
sprays  in  air  washers  should  be 
eliminated  as  these  units  have  been 
associated  with  several  outbreaks  of 
allergic  diseases  (Arnow  et  al.  1978) 
as  water/microbial  aerosols  migrate 
from  water  spray  systems  into  duct- 
work and  occupied  spaces. 

Operate  cooling  systems  to  provide 
40  to  60  percent  relative  humidity 
(RH)  all  year  if  possible.  If  neces- 
sary, reheat  with  heating  coils  even 
during  summer  to  stay  within  this 
range.  Many  pathogenic  bacteria  do 
not  like  the  50  percent  RH  and  will 
expire.  Relative  humidity  in  occu- 
pied spaces  should  certainly  not  ex- 
ceed 70  percent  to  avoid  fungal 
spore  germination  and  proliferation 
(Brundrett  and  Onions  1980).  Inap- 
propriately high  relative  humidity 
will  lead  not  only  to  condensation 
on  cool  surfaces  but  also  will  allow 
hygroscopic  materials  in  the  envi- 
ronment (e.g.,  human  skin  scales)  to 
absorb  enough  water  to  facilitate  mi- 
crobial growth.  Cooling  coils  of 
AHUs  should  be  run  at  a  low 
enough  temperature  to  dehumidify 
conditioned  air.  Stagnant  water 
should  not  be  allowed  to  accumulate 
under  cooling  coils  in  AHUs.  Proper 
inclination  and  continuous  drainage 
of  drain  pans  is  necessary. 
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Clean-up  of  existing  contamination 

Once  contamination  has  occurred,  two  ap- 
proaches are  available;  air  cleaning  and  sur- 
face and  source  decontamination.  Clean-up  of 
existing  contaminated  HVAC  systems  and  occu- 
pied spaces  usually  consists  of  inactivating  the 
air-conditioning  system  and  applying  disinfec- 
tants to  all  contaminated  building  spaces  as 
well  as  ductwork,  coils,  inside  surfaces  of 
AHUs,  fans  and  other  equipment  in  the  air 
stream.  Washing  microbially-contaminated 
hard  surfaces  with  5-10  percent  bleach  solution 
will  kill  resident  microorganisms  but  will  not 
prevent  recontamination  unless  underlying  con- 
ditions are  changed.  Other  cleaning  procedures 
include  the  use  of  biocides.  These  procedures 
often  pose  safety  questions.  Many  biocides  used 
in  indoor  environments  will  enter  the  air  of  the 
occupied  space  and  may  pose  greater  risks  than 
the  microorganisms  that  they  are  designed  to 
kill.  New  biocides  that  are  firmly  attached  to 
the  treated  substrate  demonstrate  promise  but 
have  to  be  adequately  tested  before  use  on  a 
large  scale  (Burge  1990).  Biocide  activity  in- 
creases with  temperature  and  should  be  ap- 
plied for  a  period  of  one  to  three  hours, 
depending  on  its  type.  Biocide  types  include: 
Phenolics,  alcohol,  hydrogen  peroxide,  iodine, 
glutaraldehyde,  quaternary  ammonium 
(quats),  hypochlorites  and  formaldehyde  (a  pos- 
sible carcinogen).  There  are  many  other 
biocides  which  are  used  for  special  applications 
such  as  cooling  towers  (halogenated  nitro- 
propane,  ethylene  dichlorites,  carabamates  and 
sodium  hypochlorites).  The  biocide  should  be 
registered  and  meet  FDA  and  EPA  standards. 
Biocides  are  regulated  by  EPA's  Environ- 
mental Pesticide  Act. 

Cleaning  techniques  vary  widely  but  can  be 
generalized  as  follows: 

•  Inactivate  the  HVAC  system. 

•  Vacuum  with  a  cleaner  and  equip- 
ment which  discharges  outside  the 
building. 

•  Remove  contamination  by  cleaning 
with  a  detergent. 

•  Apply  steam  cleaning  to  contami- 
nated areas. 


•  Apply  biocide  in  accordance  with  di- 
rections approved  by  the  appropri- 
ate agency.  Be  certain  that  the 
biocide  is  not  aerosolized  and 
spread  to  occupied  spaces. 

•  Flush  with  an  inactivator  after  the 
proper  application  interval  is  com- 
pleted. 

"New"  common  sense  parameters: 

Eliminate  all  underground  ductwork.  Under- 
ground ducts,  along  with  ducts  constructed  of 
organic  or  insulation  type  materials,  have  been 
the  cause  of  serious  contamination  conditions. 
Ducts  can  easily  become  contaminated  with 
pathogens  (e.g.,  hospital  operating  rooms)  in  a 
wet,  underground  environment,  and  are  practi- 
cally impossible  to  clean.  Ductwork  within 
buildings  may  be  less  prone  to  such  contamina- 
tion. If,  however,  underground  ductwork  can- 
not be  eliminated,  special  attention  should  be 
given  to  moisture  and  substrate  control  in  duct- 
work. 

Preventive  maintenance.  The  ventilation  sys- 
tem should  be  carefully  examined  for  appropri- 
ate design,  operation  and  maintenance.  For 
example,  indoor  microbial  contamination  can 
result  when,  for  energy  conservation  purposes, 
chilled  water  systems  are  run  at  temperatures 
too  high  to  provide  adequate  dehumidification. 
Likewise,  "economizer"  systems,  designed  to  in- 
crease the  amount  of  outside  air  brought  into 
the  system  during  mild  weather,  may  be  set  so 
that  more  than  minimum  fresh  air  is  delivered 
during  hot  and  cold  weather.  Keeping  all  parts 
of  the  ventilation  system  clean  is  of  paramount 
importance. 

If  the  indoor  environment  is  kept  clean  and 
ductwork  is  carefully  cleaned  by  experienced 
workers,  conditions  necessary  for  microbial 
growth  cannot  occur.  The  AHUs  should  be  con- 
structed or  modified  so  that  maintenance  per- 
sonnel have  easy  access  to  the  heat  exchange 
components  as  well  as  to  drain  pans.  Access 
doors  should  be  strategically  placed  in  all 
AHUs.  The  emphasis  should  be  placed  upon 
cleaning  ducts.  Although  many  buildings  re- 
ceive excellent  daily  housekeeping,  ducts  and 
AHUs  may  not  have  been  touched  since  the 
HVAC  equipment  was  installed.  Filters  usually 
receive  attention,  at  least  every  six  months,  but 
coils  are  often  partially  clogged,  green  or  black 
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mold  and  fungi  may  be  growing  rapaciously  on 
the  coils  or  drain  pan.  Aerosolization  can  easily 
occur  when  the  growth  is  sufficient,  and  if  the 
fungus  is  a  facultative  parasite  ( i.e.,  Aspergil- 
lus, Cryptococcus  or  Mucor),  it  can  use  any 
other  living  organism  (even  a  person)  as  a  host. 
Even  if  contamination  of  occupants  has  not  oc- 
curred, fungi  produce  metabolic  substances 
that  adversely  affect  the  indoor  environment. 
Biocidal  mycotoxins  are  carcinogenic  fungal  me- 
tabolites. In  addition,  as  indicated  earlier, 
fungi  can  produce  toxic  chemicals  and  emit 
VOCs  as  they  digest  substrates. 

The  maintenance  program  should  also  in- 
clude the  regular  checking  of  all  temperature 
(and  humidity)  control  settings  and  their  per- 
formance; routine  servicing  of  the  HVAC  plant; 
and  regulation  of  fresh  air  quantities,  air  move- 
ment and  air  distribution. 

Test  for  bioaerosols  frequently.  Such  tests  fo- 
cus on  environmental  determinants  of 
bioaerosols,  including  contributions  from  out- 
door aerosols,  ventilation  mode  and  rate,  and 
factors  allowing  intrusion  or  accumulation  of 
moisture  indoors.  The  nature  of  the  landscape 
can  affect  levels  of  airborne  spores  immediately 
surrounding  a  building,  especially  low-rise  and 
large  buildings  with  air  intakes  near  the 


ground.  Dense  vegetation  close  to  a  building 
may  increase  the  risk  of  microbial  contamina- 
tion within  (Kozak  et  al.  1978).  Also,  farming 
activities  produce  dense  bioaerosols  that  can 
penetrate  nearby  buildings  (Burge  1990). 

Install  appropriate  filters.  For  buildings  that 
are  contaminated  with  microbial-infected  mate- 
rial, it  is  better  to  discard  the  contaminated 
items  (carpet,  upholstery,  ceiling  tiles,  etc.) 
rather  than  to  attempt  to  disinfect  them;  how- 
ever, if  the  extent  of  contamination  is  limited 
or  uncertain,  disinfection  may  be  attempted  us- 
ing high-efficiency  particulate  air  (HEPA)  fil- 
ters. These  filters  are  introduced  into  the 
airstream  immediately  before  it  enters  the 
space  so  that  any  contaminants  generated 
within  the  system  itself  are  also  caught  (Figure 
5).  Since  it  is  almost  impossible  to  clean  con- 
taminated suspended  ceiling-return  air 
plenums,  these  may  have  to  be  by-passed  by  in- 
stallation of  return  air  ductwork.  During  build- 
ing clean-up,  microbially  laden  materials 
should  be  carefully  removed  so  as  to  minimize 
aerosolization  of  inhalable  particulates. 

The  selection  of  filters  varies  according  to  in- 
itial and  maintenance  costs,  energy  costs  and 
perceived  levels  of  "quality"  required  for  a  build- 
ing space.  Generally,  as  air  filtration  efficiency 
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Figure  5.  Location  of  HEPA  filter  in  an  air  handling  unit 
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Figure  6.  Efficiency  of  high-performance  dry-media  filters 


increases,  the  initial  cost  of  the  filter  increases, 
and  the  air  pressure  drop  across  the  filter  in- 
creases. An  increase  in  air  pressure  drop  re- 
sults in  an  increase  in  the  required  fan  power 
and,  thus,  an  increase  in  operating  cost.  Figure 
6  and  Table  2  show  the  efficiency  of  different 
high  performance  filters  as  function  of  particle 
size. 

One  type  of  air  cleaning  device — the  electro- 
static air  cleaner — has  both  reasonable  effi- 
ciency and  a  low  pressure  drop,  although  the 
initial  cost  is  high.  In  this  device,  air  is  drawn 
through  an  electrostatic  field  in  which  particles 
in  the  air  take  on  an  electrostatic  charge  and 
are  attracted  to  electrostatically  charged  plates 
where  they  are  collected  (Figure  7).  An  advan- 
tage of  these  types  of  air  cleaners  is  that  they 
are  effective  even  on  very  small  particles.  One 
drawback  of  these  devices,  however,  is  the  gen- 
eration of  small  quantities  of  ozone  which  may 
present  a  health  hazard.  Another  disadvantage 
is  that  as  the  build-up  of  dust  assembles  on  the 
charged  plates,  globules  of  charged  dust  some- 
times spall  off  and  continue  downstream  until 


Table  2. 

Engineering  Data  of  High  Performance  Dry- 
(corresponding  lo  Figure  6) 

Media 

Fillers 

12*24*8 

Dimension  (inch) 
12*24*12                      24*24*8 

24*24*12 

Media 
Type 

Rated  Capacity 

(ft3/min) 

Pressure  Loss 
(in.  H2O) 

M-2A 
M-15 
M-100 
M-200 

900 
900 
650 
450 

1025                                  1725 

1025                                  1725 

875                                  1325 

630                                    920 

2000 
2000 
1700 
1200 

0.15 
0.35 
0.40 
0.40 

14.5 

Effective  Filtering  Area  (ft2) 
20.8                              29.0 

41.7) 

Note:  1  in.  =  25.4  mm; 


1  in.2  =  645.2  mm2; 


lft7min  =  0.472  1/s; 


1  in.  H2O  =  248.8  Pa 
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(a)  Air  cleaner  inserted  in  duct 


Alternate 
4  j plates  grounded 


-  [^^^^J^l>^:;:^v-A/^/l■v^^^■■-^l^T>■v>J^w^^A:^^/>:^^w^^/J.r>»v'Ay.■^^^^v>/?gl    — 
\    \  I    i   / 

\\j  //Airflow  c 

— v//l\\  ^- Intermediate  plates 

7  / 1  \  \  '  *+  charged  to  high 


j!r<A2tX3*ratt.tzTm-s3X?.i?A;s)si/3*?4),r\.,';/zs\3'SJ/7Wiy<^ 


EJEB^B^SSZaaSESB^S^aggBB   -       Positive  potential 


Wires  at  high  -^feg^,, ,f^,^l  _      of  charged  dust 

//|\\  -=zZ^~-^-^- particles 


-i-    Theoretical  paths 


positive  potential 


//!\\    /' 


»r^.»j«m*^^v^^VA^ibv?^.w«xr^iiM«Mi»vMrj-riicjx.rj».iir.s^a-CTv«r^i.'- 


Positively  charged 
particle 


^ 


(b)  Detail  showing  components  and  charging  mechanism 


Figure  7.  Electrostatic  air  cleaner 
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caught  by  the  cooling  coils  or  the  dust  is  distrib- 
uted to  the  occupied  space.  If  this  happens  on  a 
large  scale,  a  considerable  build-up  of  space 
charge  will  result,  which  not  only  drives  the 
charged  dust  particles  to  the  wall  and  causes 
blackening  of  the  surface,  but  it  may  represent 
a  certain  health  hazard  to  space  occupants. 
Electrostatic  fdters  can  be  automatically 
washed  to  eliminate  spall-off,  but  there  is  no  as- 
surance that  spalling  will  not  occur  in  between 
washing  cycles. 

Instead  of  integrating  expensive  high-effi- 
ciency fdtering  systems  with  the  central  HVAC 
unit,  another  alternative  would  be  the  zoning 
of  stand-alone  air  fdtering  equipment  operat- 
ing independently  of  the  central  HVAC  system 
to  achieve  both  a  high  circulation  rate  and 
proper  air  mixing.  No  ductwork  is  required. 
Thus,  each  zone's  filtering  system  has  its  own 
fan,  which  can  operate  either  with  or  without 
the  central  HVAC  fan.  Air  circulation  through 
these  filters  should  occur  at  rates  between  6 
and  10  space  air  changes  per  hour. 

Design  properly  maintainable  mechanical  sys- 
tems and  spaces.  A  frequent  cause  of  microbial 
contamination  is  poor  engineering  design  of  the 
system  leading  to  difficulties  in  getting  access 
to  certain  components  of  the  HVAC  system. 
The  following  examples  illustrate  such  cases: 

•  AHUs  are  designed  without  access 
doors  to  the  heat  exchanger  section. 

•  Small  AHUs  and  heat  pumps  are  lo- 
cated in  inaccessible  spaces  above 
ceiling  tiles  prohibiting  entry  into 
mechanical  components  to  clean 
coils,  drain  pans,  etc. 

•  Some  mechanical  rooms  and 
plenums  are  so  confining  that  ac- 
cess is  rendered  difficult  if  not  im- 
possible once  the  equipment  is 
installed.  Architects  must  be  made 
to  provide  adequate  space  when 
buildings  are  designed  or  remodeled. 

In  addition  to  remedying  the  above,  fan  coil 
unit  designs  should  facilitate  exposure  of  drain 
pans,  cods,  etc.,  as  quickly  as  possible.  Figure  8 
illustrates  the  configuration  of  equipment  of  a 
generic  HVAC  system  for  a  building  experienc- 
ing VOCs,  airborne  infectious  diseases  and 


chemical  and  particulate  contamination.  The 
"worst-case  scenario"  is  not  likely  to  occur  ex- 
cept in  severe  cases  of  contamination.  A  high 
level  of  radiation  can  be  emitted  from  ultra-vio- 
let (UV)  lamps,  located  in  the  AHU  adjacent  to 
the  carbon  canisters,  and  it  is  vitally  important 
to  interlock  the  lamps  with  the  access  door 
latch  and  the  fan  motor  to  prevent  eye  injury  to 
operating  engineers.  While  the  cost  of  such  a 
system  may  be  high,  it  will  most  likely  be  lower 
than  the  cost  of  decontaminating  an  infected 
building.  Moreover,  for  new  construction,  is  it 
prudent  to  invest  in  an  on-line  decontamina- 
tion system  up  front  or  take  the  chance  that 
the  building  will  escape  the  fate  of  many?  Once 
a  building  is  contaminated,  remedial  actions 
can  be  extremely  expensive  and  the  loss  of  ten- 
ants (as  well  as  reputation)  is  sometimes  unre- 
coverable. 

CONCLUSION 

The  major  contamination  experienced  by 
most  HVAC  systems  is  microbial  in  the  form  of 
a  bioaerosol.  Since  most  bacteria  and  fungi  are 
facultative  parasites  and  grow  on  both  living 
and  nonliving  organic  material,  it  is  possible 
for  an  air-conditioning  system  to  act  as  an  am- 
plifier or  disseminator.  This  applies  especially 
to  saprophytic  microorganisms,  fungi,  protozoa 
and  most  bacteria.  With  the  exception  of  Le- 
gionella which  must  reside  in  a  water  droplet, 
if  the  substrate  hosting  the  organism  is  dis- 
turbed, the  organism  becomes  airborne.  Thus, 
the  air-conditioning  system  becomes  an  accom- 
plice in  bacterial  dissemination. 

While  we  have  little  control  of  the  external 
environment,  we  can  do  something  about  the  in- 
door environment.  Remediation  can  be  in  the 
form  of  a  clean-up  (after  contamination  is  dis- 
covered) or  on-line  which  assigns  remediation 
to  the  air-conditioning  system  while  it  is  operat- 
ing. A  manual  clean-up  would  mean  the  air-con- 
ditioning system  would  be  deactivated  for  a 
period  of  time  while  the  cleaning  and  disinfec- 
tion is  done.  Should  the  role  of  the  air-condi- 
tioning system  be  expanded  to  include  air 
quality  control?  If  so,  definition  of  the  new  de- 
sign parameters  for  such  air-conditioning  sys- 
tems should  be  assigned  and  adopted  by 
ACGIH,  ASHRAE,  BOCA  and  others,  as  recom- 
mendations, standards  or  code  requirements. 
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The  following  list  specifies  equipment  and 
performance  of  a  generic  HVAC/IAQ  system  to 
detect  and  remove  bacterial  contamination  of 
all  types.  Design  parameters  for  the  new  or  ret- 
rofit, on-line  air  quality  control  system  would 
include: 

•  Conformance  during  the  commis- 
sioning process  with  recommended 
microbial  exposure  range  limits. 
Based  on  the  cited  literature,  the  up- 
per limit  should  not  exceed  1000 
CFU/m3.  Reference  is  also  made  to 
the  ACGIH  (1989)  guidelines  recom- 
mending comparison  of  indoor  and 
outdoor  concentrations. 

•  Airborne  hypersensitivity  disease 
detection. 

•  Biocidal  treatment  of  the  air  stream 
while  diverting  the  discharge  air 
outside  the  building,  and  expiration 
by  UV  lamps  of  certain  other  active 
bioaerosols. 

•  Depyrogenation  cycles  through  the 
AHU  at  540°F  (280°C)  while  divert- 
ing the  air  outside  (filters  should  be 
removed  during  this  cycle  to  pre- 
vent damage). 

•  New  standards  for  equipment  manu- 
facture to  include  materials  and  as- 
sembly. AHUs  and  ductwork  should 
not  have  insulation  on  the  inside 
surface  which  can  collect  and  am- 
plify microorgansims;  inside  duct 
surface  should  be  smooth  such  as 
aluminum  or  stainless  steel;  water 
atomization  humidifiers  should  not 
be  used  (steam  only);  all  equipment 
designed  for  high  temperature 
depyrogenation  if  this  cycle  is  in- 
cluded; ducts  should  have  flanged 
joints;  gaskets  and  sealants  should 
not  contain  organic  or  VOC  pollut- 
ants. 

•  A  microbial  sensing  system  can  pos- 
sibly be  installed  to  detect  a  build- 
up in  concentration  passing  through 
the  AHU.  The  bacterial  colony  mor- 
phology should  be  established  at  ge- 
neric level  {Bacillus,  Strepococcus, 


Micrococcus  and  Corynebacterium) 
for  comparison  with  test  samples. 

•  As  previously  stated,  ductwork 
should  not  be  installed  below 
ground,  if  possible.  Existing  under- 
ground ducts  should  be  removed  or 
filled  with  a  strong  biocide  and  sand 
combination  and  sealed  with  a  thick 
concrete  topping  layer,  caulked  and 
sealed  at  the  edges.  New  return 
ducts  should  be  installed  for  exist- 
ing systems.  (If  ceiling  diffusers  are 
used,  use  low  sidewall  return  regis- 
ters). 

•  Buildings  should  have  positive  pres- 
sure with  respect  to  the  outside. 

•  A  bake-out  cycle  could  possibly  be 
designed  to  provide  100  to  105°F 
(38  to  41°C)  air  throughout  the 
building  to  drive  off  VOCs  and  re- 
duce bioaerosols.  The  bake-out 
equipment  should  include  two  run- 
around  coils,  piping  and  a  small 
pump  to  reduce  operating  costs.  The 
initial  bake-out  should  occur  after 
the  building  is  completed  and  before 
occupancy.  Only  100  percent  out- 
side air  should  be  used,  and  this  air 
should  be  entirely  exhausted.  The 
high  temperature  air  should  not  be 
continuously  circulated  throughout 
the  building,  but  should  take  advan- 
tage of  thermal  "pumping"  charac- 
teristics to  induce  VOC  release  from 
fungi  and  building  materials,  i.e.,  2 
hours  at  105°F  (41°C),  3  hours  at 
80°F  (27°C),  repeated  over  a  48 
hour  period.  Continuous  105°F 
(41°C)  air  within  a  building  for  two 
days  is  likely  to  cause  extensive 
damage  to  the  building  such  as 
warping  of  doors,  wrinkling  of  car- 
peting, possible  fracture  of  windows 
and  thermal  expansion  of  various 
components.  The  building  should 
also  be  brought  up  to  temperature 
slowly— by  10°F  (6°C)  increments 
over  a  12  hour  preheat  period  [10°F 
(6°C)  increase  every  4  hours]. 

Other  air-conditioning/quality  system  re- 
quirements would,  of  a  necessity,  include  re- 
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moval  of  ozone,  VOCs,  chemicals  and  particu- 
lates, subjects  not  part  of  this  discussion.  While 
it  is  unlikely  that  all  nine  requirements  would 
be  necessary  for  any  one  building,  the  above  pa- 
rameters can  be  applied  individually  as  re- 
quired. 

It  appears  that  the  building  air-conditioning 
system,  which  has  remained  basically  un- 
changed through  the  years,  may  become  a 
"building  air  quality  control  system."  Much  re- 
search is  necessary  to  determine  how  we  can 
use  building  HVAC  systems  to  control  the  in- 
door environment,  preserve  our  health,  and  pro- 
tect us  from  sickness  generated  in  the 
workplace.  HVAC  systems  must  function  to  pro- 
vide an  integrated  environment  where  tempera- 
ture, air  distribution,  humidity,  energy 
consumption  and  IAQ  are  maintained  within 
acceptable  levels. 
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ABSTRACT 


The  complexity  of  comprehensive  repre- 
sentations for  building  performance  simulation 
is  illustrated,  based  on  the  application  of 
META-4  (mass  and  energy  transfer  analysis)  in 
multi-layered  building  enclosure  components. 
Required  multiple  modeling  strategies  (nu- 
meric computation,  geometric  procedures,  rule- 
based  approaches,  generative  and  "two-way" 
systems)  are  discussed.  Frameworks  for  effec- 
tive design  supporting  communicative  struc- 
tures are  proposed. 
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INTRODUCTION 


Continuous  advancements  in  computer  tech- 
nologies have  spurred  the  development  of  nu- 
merous computer  aided  design  tools  to  assist  in 
performance-based  design  of  buildings.  This 
generally  involves  the  analysis  and  evaluation 
of  those  design  features  that  affect  the  function- 
ality and  occupational  quality  of  a  building. 
However,  many  of  these  tools  have  found  lim- 
ited application  in  architectural  practice.  Even 
then,  normally  only  specialist  professionals  use 
these  performance  simulation  models,  often  for 
verification  purposes;  there  appear  to  be  three 
main  reasons  for  this. 

First,  the  tools  have,  by  and  large,  been  de- 
veloped as  stand-alone  applications  within  a 
specific  design  domain.  The  underlying  pro- 
gram structure  (e.g.,  process  and  data  manage- 
ment) reflects  certain  idiosyncratic  character- 
istics that  make  it  difficult  for  integration  with 
other  essential  performance  evaluation  and  geo- 
metric modeling  tools  (Mahdavi  and  Lam 
1991).  Second,  the  problem  of  effective  informa- 
tion transfer  between  the  user  and  the  compu- 
tational environment  is  still  unresolved.  Third, 
efforts  in  tool  development  conceptually  fall 
short  of  challenging  the  existing  (generally  se- 
quential and  rather  fragmented)  building  deliv- 
ery process  (Mahdavi  and  Lam  1993). 

This  paper  aims  to  highlight  the  complex 
representations  involved  in  building  perform- 
ance simulation  and  to  discuss  various  model- 
ing strategies  adopted  in  the  META-4  (Mass 
and  Energy  Transfer  Analysis  for  multi-lay- 
ered building  enclosure  components)  computa- 
tional environment,  which  include  numerical 
computation,  geometric  procedures,  rule-based 
approaches,  as  well  as  generative  and  "two- 
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way"  systems.  Suggestions  for  appropriate  com- 
municative frameworks  are  offered. 

CRITIQUE  OF  PARTIAL  SPECIFIERS 

In  the  past,  the  limitations  of  tools  in  terms 
of  computational  capacity  and  speed  made  it 
necessary  to  reduce  the  quantity  of  information 
processing  for  design  decision  making  by  intro- 
ducing "partial  specifiers"  and  rules  of  thumb. 
Given  the  increased  computational  capabilities 
in  building  performance  simulation,  many  of 
these  specifiers  and  rules  should  be  critically 
reviewed.  For  instance,  simplified  performance 
indicators  such  as  U-values  and  Sabine  rever- 
beration times  are  used  as  exemplification  of 
performance  representations  ("structures  for 
deriving  thermal  and  acoustical  performance") 
in  deriving  object-oriented  engineering  design 
databases  (Eastman  et  al.  1991). 

Aside  from  the  fact  that  U-value  modeling 
cannot  represent  transient  heat  transfer  phe- 
nomenon, various  research  studies  (Mahdavi 
1993)(Mahdavi  et  al.  1992)(Panzhauser  et  al. 
1990)  have  demonstrated  the  limitations  of  U- 
value-based  calculations  for  the  determination 
of  heat  flows,  surface  temperatures,  and  criti- 
cal relative  humidities  (essential  for  the  evalu- 
ation of  surface  condensation  risk)  under 
steady-state  conditions.  The  underlying  as- 


sumption of  U-Value-based  calculations  is  in- 
sufficient to  deal  with  many  of  the  complex  com- 
posite constructions  (which  involve  thermal 
bridging)  and  could  yield  erroneous  results  in 
both  transmission  heat  loss  and  interior  sur- 
face temperature  calculations. 

Figure  1  illustrates  the  point  by  comparing 
simulation  results  of  minimum  interior  surface 
temperatures  using  both  the  detailed  numeric 
method  (based  on  finite  difference  computation 
of  conductive  heat  transfer)  versus  the  simpli- 
fied U-value-based  method  for  30  planar  build- 
ing components  typical  of  residential  construc- 
tions in  North  America  (Mahdavi  and  Mathew 
1992).  In  this  case,  the  simplified  method  tends 
to  provide  over-estimated  values  in  most  in- 
stances. 

In  acoustics,  notwithstanding  the  inherent 
shortcoming  of  the  "classical"  notion  of  rever- 
beration time,  the  empirical  Sabine  algorithm 
for  computing  reverberation  time  is  also  sub- 
ject to  limitations  as  it  is  valid  only  for  rooms 
with  low  average  absorption  (a,,,  <  0.2)  and,  con- 
trary to  expectation,  it  gives  reverberation 
times  above  zero  for  o^  =  1  (ep.  Figure  2). 

Similarly,  the  Sound  Reduction  Index  (Rw)  of 
the  building  elements  (walls,  floors)  has  tradi- 
tionally been  used  to  specify  sound  transmis- 


10  12  14  16  18  20 

ts„min  [°C]      (Numeric  Method) 


Figure  1 .  Comparative  analysis  of  indoor  surface  temperatures 
for  30  residential  constructions  using  different  simulation  meth- 
ods. (Assumed  boundary  conditions:  Outdoor  temperature  =  - 
1  5°C,  Indoor  temperature  =  20°C.) 


Figure  2.  Comparative  illustration  of  reverberation  time  compu- 
lation according  to  Sabine  and  Eyring  for  3  different  room  con- 
figurations with  S=  total  room  surface  area,  V=  room  volume. 
[1:  S=130  m2  (1  400  ft2),  V=100  m3  (3532  ft3);  2:  S=600  m2 
(6458  ft2),  V=  1000  m3  (3 53 15  ft3);  3:  S=2800  m2  (30139 
ft2),  V=  10000  m3  (3531  47  ft3).] 
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sion  between  adjacent  spaces.  However,  the 
sound  transmission  is  not  only  a  function  of  the 
sound  insulation  of  the  partition  element  but 
also  the  result  of  the  combined  effects  of  sound 
transmission  through  flanking  paths  (ep.  Figure 
3).  Recent  documents  have  identified  the  nor- 
malized sound  level  difference  (Dn_T|W)  as  a 
more  appropriate  indicator  for  the  formulation 
of  building  acoustic  requirements  (Mahdavi 
1991). 

A  similar  configuration  to  that  shown  in 
Figure  3  is  used  to  demonstrate  the  potential  im- 
pact of  the  flanking  transmission  on  the  overall 
sound  level  differences  (Mahdavi  and  Lam 
1991).  Figure  4  illustrates  the  deviation  of  the 


sound  reduction  index  of  the  partition  element 
from  the  normalized  sound  level  difference  be- 
tween two  adjacent  rooms.  The  relevant  under- 
lying assumptions  are  summarized  in  Table  1 
below. 

The  examples  point  to  the  need  to  recognize 
and  understand  the  complexities  involved  in 
performance  evaluation  and  to  derive  appropri- 
ate frameworks  for  accommodating  and  repre- 
senting these  complexities  in  computer-aided 
design  environments.  While  it  is  arguable  that 
particularly  in  the  early  design  stage  a  higher 
level  of  abstraction  is  acceptable,  it  does  not 
mean  that  rigor  and  accuracy  in  the  simulation 
process  can  be  compromised. 


Table  1 

.  Matrix  of  modeling  assumptions 

(cp. Figure  4) 

Element 

Area 

[m2] 

Surface  Density 
[kgm2] 

Rw 

[dB] 

Junction 
Geometry  [+,  T] 

Partition  Wall 

12 

100  to  700 

38.8  to  66.2 

N/A 

Exterior  Wall 

15 

400 

56.3 

T 

Interior  Wall 

15 

100 

36.8 

+ 

Floor 

20 

300 

52.3 

+ 

Ceiling 

20 

300 

52.3 

+ 

SOURCE 
ROOM 


RECEIVER 
ROOM 


70 


60 


5    55     - 


50 


o 

•i) 


40 

35  1 — i — i — i — i — i — r 

100  200  400 

Surface  Density  (kg  m': 


soo 


Figure  3.  Direct  and  flanking  sound  propagation  paths  be- 
tween adjacent  spaces 


Figure  4.  Deviation  of  the  sound  level  difference  between  two 
adjacent  rooms  in  a  typical  massive  construction  from  the 
sound  reduction  index  of  the  partition  element  as  a  function  of 
the  surface  density  of  the  partition  element  (cp.  Table  1  for 
modeling  assumptions) 
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MULTIPLE  STRATEGIES  FOR  COMPLEX 
REPRESENTATIONS 

To  illustrate  the  complexity  involved,  the  de- 
velopment of  a  computational  module  for  simu- 
lating and  visualizing  heat  and  moisture 
transfer  (META-4)  in  a  multi-layered  building 
enclosure  component  is  presented.  This  pro- 
gram mainly  addresses  the  need  for  up-front 
evaluation  of  building  enclosure  components, 
particularly  with  respect  to  their  long-term  hy- 
gro-thermal  behavior  (Mahdavi  and  Lam  1993). 

As  the  following  sections  demonstrate,  a  rig- 
orous computation  of  performance-related  phe- 
nomena may  include  comprehensive 
implementation  of  multiple  strategies  whereby 
a  variety  of  approaches  such  as  numeric  compu- 
tation, geometric  procedures,  rule-based  sys- 
tems, generative  and  "two-way"  routines 
(Mahdavi  and  Berberidou-Kallivoka  1992)  as 
well  as  data  management  techniques  and 
graphic  information  representation  are  utilized 
simultaneously. 

Application  of  Numeric  Methods 

Treatment  of  Boundary  Conditions.  In  estab- 
lishing the  external  boundary  conditions,  it  has 
been  common  practice  to  adopt  the  dry  bulb  air 
temperature  from  the  available  hourly  weather 
data  files.  This  does  not  take  into  account  the 
significant  impact  of  solar  radiation  on  surface 
temperatures  of  the  building  envelope. 
Weather  files  generally  do  not  include  solar  ir- 
radiance  on  arbitrarily  oriented  surfaces.  To  ob- 
tain these  data  and  to  estimate  the  intensities 
of  direct  and  diffuse  solar  radiation  on  an 
hourly  basis,  a  solar  radiation  pre-processor 
has  been  incorporated  in  META-4  (Heindl  and 
Koch  1976). 

Based  on  the  radiation  intensities,  the  "ra- 
diation air  temperature"  is  then  modeled  ac- 
cording to  Equation  1  (Koch  and  Pechinger 
1977),  which  can  be  solved  utilizing  the  Newto- 
nian Approximation  Theorem. 

ak(TL-T')  +  (as-J)  +  5.18- 

lO-8^  [euEuT4u  +  (l-eu)egT4L]- 

[5.78    10-8-eo(T*)4]  =  0  (1) 


a^  =  convective  heat  transmission  coefficient  in 
W-m^-K1  (Btu-h^-ff  2-R_1) 

Tj,  =  outdoor  air  temperature  in  K  (R) 

T   =  radiation  air  temperature  in  K  (R) 

as  =  absorption  coefficient  of  building  surface  for 
solar  radiation 

J  =  solar  radiation  intensity  in 
W-m-2  (Btu-h_1-ft"2) 

e0  =  emittance  of  building  surface 

eu  =  angle  factor  of  surrounding 

eu  =  emittance  of  surrounding 

eg  =  emittance  of  atmosphere 

Tu  =  surrounding  temperature  in  K  (R). 

It  is  also  conceivable  that  internal  boundary 
conditions  will  also  fluctuate  depending  on  occu- 
pancy and  mechanical  system  schedules  and 
settings.  A  procedure  has  also  been  incorpo- 
rated in  the  program  to  deal  with  this  scenario. 

FD-based  Heat  Transfer  Modeling.  The  prin- 
ciple governing  the  heat  conduction  equation  is 
based  on  Fourier's  Law: 


Qx  -  A  (7»  '  (  -   g  x  ) 


(2) 


where 


q^  =  heat  transfer  rate  in  the  x  direction  in 


W-m"2  (Btu-h-1-ft-2) 


•1  xr-l 


\(Tu)  =  thermal  conductivity  in  W-m'   K 

(Btu-h'^ft^R"1)  f°r  temperature  T  in  K  (R) 
and  moisture  content  u  (%). 

Under  steady-state  conditions  (with  tempera- 
ture difference  AT),  the  temperature  distribu- 
tion within  a  multi-layered  planar  building 
element  (with  n  layers  of  thickness  dj  and  con- 
ductivity Xj )  can  be  derived  from: 


where 


Q  =  (R»  +  X  djXj  + 

V  7=1 


R* 


■AT 


(3) 
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where 

q  =  heat  transfer  rate  in  W-m"2  (Btu-h-1-ft"2) 

2  1 

Rs:-  =  inside  surface  resistance  in  m  -K-W" 
(ft"2-R-h-Btu-l) 

Rsc  =  outside  surface  resistance  in  m  -K-W" 
(ft^-R-h-Btu"1). 

As  with  the  limitations  of  the  U-value  speci- 
fier, the  steady-state  simulation  procedure  does 
not  address  the  dynamic  nature  of  the  bound- 
ary conditions  and  its  impact  on  the  heat  trans- 
fer rate.  Transient  simulation  procedures  are 
now  regarded  as  essential  for  any  thermal 
evaluation.  To  achieve  a  higher  level  of  resolu- 
tion in  the  hygro-thermal  analysis,  each  mate- 
rial layer  of  the  building  component  is  initially 
discretized  into  a  predefined  number  of  sub- 
layers. An  additional  discretization  criterion 
limits  the  maximum  thermal  resistance  value 
of  each  sub-layer  so  as  to  maintain  the  desired 
computational  resolution. 

The  "effective"  surface  resistance  Rs*  is  cal- 
culated as  a  function  of  the  radiative  heat 
transfer  coefficient  a.  and  the  convective  heat 
transfer  coefficient  ak: 


Use    =  (as  +  a-k ) 


[m2-K-W  l] 


(4) 


Referring  to  Equation  1,  taking  the  deriva- 
tive relative  to  71*,  the  radiative  heat  transfer 
coefficient  can  be  determined  as  follows: 

a*=F(T*)=  -<xk- 


[5.78  ■  If/8  •  e0  •  4  (  T  ) 3  ] 


[W-m  2K  l  ] 


(5) 


The  convective  heat  transfer  coefficient  ak  is 
derived  as  a  function  of  the  wind  velocity  v 
(m-s  ): 


a*  =  5.9  +  4.1  •  v  [W-m  2-K  1  ] 


(6) 


The  "effective"  L/*-value  of  the  enclosure  compo- 
nent is  then  given  by: 


tf*  =  (««  +  £  djXj  +  Rse')~l 


;=i 


[W-m~2-K  *] 


(7) 


The  calculation  of  transient  temperature  dis- 
tribution in  the  sub-layers  is  based  on  the  sub- 
stitution of  the  transport  equation  by  finite 
differences  and  the  corresponding  balance  equa- 
tions are  substituted  by  finite  control  volumes 
and  solved  by  the  tri-diagonal  matrix  solution 
technique.  Latent  heat  exchange  associated 
with  condensation  of  water  vapor  or  evapora- 
tion of  moisture  in  the  component  sublayers  is 
accounted  for  in  the  procedure  and  expressed 
as  a  heat  source  or  heat  sink  term  in  the  finite 
control  volume  formulation. 

Dynamic  Modification  of  Thermal  Conductiv- 
ity. A  procedure  has  been  implemented  to  dy- 
namically model  the  effect  of  moisture  content 
on  the  thermal  conductivity  of  the  materials 
(Sigmund  1984)(Wiese  1985).  The  thermal  con- 
ductivity (A.)  for  the  material  with  a  given  mois- 
ture content  (u)  can  be  derived  from  the 
thermal  conductivity  of  the  material  {Xde/)  with 
the  "default"  moisture  content  (u^): 


X=%def  ■(l  +  Z)-(l+Zdef  f1 

[W-m^-K1  ] 


(8) 


In  this  equation,  Z  is  the  thermal  conductiv- 
ity incremental  factor  related  to  the  moisture 
content  of  the  material.  Z^  refers  to  the  incre- 
mental factor  at  default  moisture  content, 
which  is  a  statistically  defined  value  below 
which  90%  of  the  empirical  measurement  re- 
sults lie  (Gosele  and  Schiile  1983).  Volume  and 
mass-related  default  moisture  content  for  vari- 
ous common  building  materials  can  be  obtained 
from  DIN  1981.  Examples  of  values  for  Z^f  are 
given  in  ONORM  1989. 

Computation  of  Water  Vapor  Diffusion.  For 
the  calculation  of  water  vapor  diffusion,  the  ma- 
terial property  u  (water  vapor  diffusion  resis- 
tance index)  of  the  component  layers  must  be 
known  (CEN  1991).  As  a  dimensionless  speci- 
fier, u  indicates  how  much  greater  the  water  va- 
por diffusion  resistance  of  the  material  is  than 
that  of  an  equally  thick  layer  of  stationary  air 
at  the  same  temperature.  Given  the  u  values, 
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Table  2.  Procedure  for  calculating  angle  a 


Condition 


a  = 


Psj>P 
Ps=P 
PsJ<P 


arctan  [}T  uy-  dj/  (psj-p  )] 
J  =  l 

90 

m 
180  -  arctan  [  £  Uj r  •  d//  (p-psj) 


the  "actual"  water  vapor  pressure  distribution 
can  be  obtained  using  the  following  equation: 


V 


(9) 


The  overall  vapor  pressure  gradient  Ap  can 
be  calculated  using  relative  humidities  {RHt , 
RHe)  and  saturation  pressures  inside  and  out- 
side (pSii ,  pSie): 


Ap  =  (RHrpsJ)-(RHeps,e) 


(ID 


where 


2       1 

g  =  water  vapor  diffusion  density  in  kg-m"  -h" 
(lb-ft"2^"1). 

Ap  =  total  vapor  pressure  differential  between 
outdoor  and  indoor  in  Pa  (psi) 

U,-  =  diffusion  resistance  coefficient  of  sub-layer  j 

(17|3)e  ,  (1/p);  =  external  and  internal  surface  dif- 
fusion resistance  in  m  -h-Pa-kg" 
(ft-2-h-psMb"1). 

For  the  internal  surface  diffusion  resistance 
(lyp); ,  a  constant  value  8300  m^h-Pa-kg"1 
(5.403  ft2-h-psi-lb"1)  for  air  temperatures  be- 
tween 10  and  20°C  (50  and  68°F)  and  tempera- 
ture difference  between  air  and  wall  of  5  to  10 
K  (9  to  18R)  is  used.  The  values  for  the  exter- 
nal surface  diffusion  resistance  (l/p)e  are  de- 
pendent on  wind  velocity.  The  following 
approximation  for  the  temperature  range  of  -20 
..  30°C  (-4  ..  86°F)  and  wind  velocity  (v)  range 
of  1  ..  20  m-s"1  (197  ..  3937  ft-min"1)  is  imple- 
mented: 


( |  \  =  100  [36  -  (120  +  60v  -  v2 )  °'5  ] 

[rn-hPakg~l]  (10) 


Application  of  Geometric  Procedures  to  Determine 
Condensation  Layer  /  Zone 

To  achieve  a  dynamic  visualization  of  the  dif- 
fusion process,  it  is  necessary  to  obtain  rapidly 
the  "true"  curve  of  water  vapor  pressure.  This 
curve  is  determined  by  adopting  a  "tangent" 
procedure  (Mahdavi  and  Lam  1993)  formulated 
as  a  generalized  solution  for  the  graphic  Glaser 
method  (Glaser  1959).  This  procedure  employs 
a  geometric  method  for  detecting  condensation 
layers  and  zones  within  the  building  compo- 
nent. Given  the  water  vapor  pressure  p  and 
saturation  curve,  an  angle  a  is  calculated  ac- 
cording to  Table  2.  Figure  5  shows  a  magnified 
portion  of  an  exemplary  building  component  to 
illustrate  this  procedure.  The  angle  a  is  calcu- 
lated sequentially  at  each  time  step  for  all  com- 
ponent sub-layers. 

Any  reversal  in  the  trend  of  a  (from  ascend- 
ing to  descending  or  vice  versa)  during  this 
mapping  process  indicates  the  start  or  end  of  a 
condensation  zone.  A  recursive  application  of 
this  tangent  procedure  identifies  the  extent  of 
condensation  layers  and  zones.  Within  the  con- 
densation zone  the  "true"  vapor  pressure  val- 
ues are  assumed  to  be  identical  to  the 
corresponding  saturation  pressures. 
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COMPONENT  SUBLAYERS 


"THEORETICAL- 
VAPOR  PRESSURE 


SATURATION  PRESSURE 


Figure  5.  Tangent  method  for  determining  "true"  vapor  diffu- 
sion path  and  condensation  layers  and  zones 


Application  of  Heuristic  Rules  to  Moisture 
Distribution 

Having  identified  the  condensation  layer(s) 
or  zone(s),  the  condensate  amount  can  be  com- 
puted. In  the  case  of  the  "block"  procedure,  long- 
term  steady-state  periods  are  defined  for  the 
process  of  condensation  and  drying.  A  dynamic 
simulation  approach  is  likely  to  provide  a  more 
accurate  representation  of  the  fluctuation  in 
moisture  content  on  a  shorter  time  scale.  The 
quantity  of  condensation  (gc)  at  a  certain  sub- 
layer m  of  a  building  component  with  n  sub-lay- 
ers is  calculated  using  the  following  equation: 

m 
go  =  (Pi  -Prn)  ■  (  1-5  •  106  •  £    Mjdj)'1    - 


(Pm-Pe)'  (1-5  -10    •£      Vjdjf 

j  =  m  +  1 


(12) 


where 


2.U-1 


gc  =  quantity  of  condensation  in  kg-m'  -h 
(lb-ft"2-*!"1) 

p;  =  internal  water  vapor  pressure  in  Pa  (psi) 

pe  =  external  water  vapor  pressure  in  Pa  (psi) 

pm  =  water  vapor  (saturation)  pressure  at  layer 
m  in  Pa  (psi). 

The  amount  of  water  igd)  drying  out  of  the 
condensation  layer  m  is  given  by: 


gd  =  (Ps-Pt)  •  ( i-s  •  io6  •  X  M/) 1  + 

n 

(ps-pe)-(1.5    106-£     u/// 

j  =  m  +  1 


(13) 


where 


g^  =  amount  of  water  drying  out  of  condensation 
layer  in  kg-m"2-h_1  (lb-ft^-h"1) 

ps  =  water  vapor  saturation  pressure  at  the  con- 
densation layer  in  Pa  (psi). 

A  procedure  currently  implemented  in 
META-4  is  based  on  "heuristic"  description  of 
the  change  in  the  mass  transfer  mechanism  as 
a  result  of  phase  change  and  the  effective  distri- 
bution mode  within  the  sub-layers  of  the  compo- 
nent. It  is  adapted  from  the  theoretical 
assumptions  used  in  the  GLASTA  program 
(Standaert  1988).  However,  the  discretization 
of  material  layers  in  META-4  provides  an  im- 
proved resolution  of  the  simulation  results.  The 
rule-based  prediction  of  the  extent  of  the  hu- 
midification  of  the  sub-layers  (L,)  is  based  on 
the  definition  of  the  critical  water  content  of  a 
material  (pWtC).  The  process  is  then  controlled 
according  to  the  following  set  of  rules: 

1-  if(Pw,i<  Pw.c.i) 

2.  then  water  vapor  diffusion  computation; 

3.  else  ifLi+1  capillary 

4.  then  transport  (pWiCii  -  pWit)  to  Li+1 

("cold  side"  sub-layer] 

5.  else  [Li+1  =  non-capillary) 

6.  transport  (pw,d  -  Pw.i) to  Wi 

("warm  side"  sub-layer} 

Statement  2  refers  to  the  mass  transfer  by 
diffusion.  Statements  3  to  6  refer  to  mass  trans- 
fer in  the  liquid  phase.  Similar  rules  are  imple- 
mented to  simulate  the  drying-out  cycle. 

Generative  Routines 

Although  the  existing  structure  of  META-4 
follows  the  conventional  "one-way"  sequence  of 
design  to  performance  transformation,  it  has  a 
clear  potential  for  adaptation  of  "two-way"  feed- 
back procedures  (Mahdavi  1993)(Mahdavi  and 
Berberidou-Kallivoka  1992).  To  demonstrate 
this  point,  the  simple  (steady-state)  condensa- 
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Diagnostics 

T 

Condensation 

Pc  >  Ps.c 

» 

Final  Layer 
Configuration 

L 

Update  Thickness 
dc  =  dc  +  Ad 

l 

T 

Update  R„  tc,  Psc 

Update  Rwv,  Pc 

Figure  6.  Simple  typical  condensation  case  within  a  multi-lay- 
ered building  component 


Figure  7.  Diagnostic  process  of  generative  modification  of 
layer  structure 


tion  case  within  a  multi-layered  building  com- 
ponent is  considered  here  (tp.  Figure  6). 

It  is  desirable  to  accommodate  routines  that 
are  activated  after  diagnostics  (e.g.  detection  of 
condensation  by  the  tangent  procedure)  to  gen- 
erate alternative  designs  that  would  avoid  con- 
densation. A  related  strategy  is  described 
below  for  the  simple  case  illustrated  in  Figure  6. 

Per  definition,  condensation  occurs  at  layer 
interface  c  if  the  computed  theoretical  vapor 
pressure  pc  is  higher  than  the  respective  satura- 
tion pressure  ps  c ,  that  is, 


Pc>Psfi 


(14) 


An  increase  in  the  thickness  of  the  layer  be- 
fore the  condensation  interface  dc  has  a  two- 
fold effect  on  the  reduction  of  the  condensation 
probability.  The  additional  thickness  and  the 
corresponding  increase  in  the  thermal  resis- 
tance Rt  leads  to  a  higher  temperature  tc  at  in- 
terface c  and  thus  to  a  higher  saturation 
pressure psc  since 


Ps.c=f(tc) 


(15) 


The  additional  thickness  also  increases  the 
vapor  diffusion  resistance  (Rwu  -  \\d),  which  in 
turn  results  in  a  lower  Pc  value: 

Pc=Pi~\-Y,  Mrf  (i  ->  c)  ■  (  X  V*h  (l  -.  n)  ] 

(16) 


Figure  7  illustrates  this  process  of  generative 
modification  of  the  layer  structure.  Of  course, 
this  procedure  must  be  extended  and  modified 
to  control  the  logical  consistency  of  the  result- 
ing thickness  and  accommodate  the  recursive 
treatment  of  the  thickness  of  other  layers  if  nec- 
essary. The  modification  should  also  address 
the  problem  of  introducing  new  layers  (e.g.,  va- 
por retarders),  multiple  condensation  zones 
and  dynamic  simulation  mode.  Despite  all 
these  limitations,  this  process  exemplifies  the 
strategic  potential  of  generative  and  "two-way" 
approaches  toward  enhancement  of  the  design 
support  capabilities  of  META-4. 

COMMUNICATION  FRAMEWORKS 

The  computational  features  of  META-4  de- 
scribed above  provide  the  potential  for  perform- 
ing complex  modeling  and  evaluation 
operations.  However,  this  potential  must  be  ef- 
fectively communicated  to  the  user  if  a  more 
pervasive  use  of  the  tool  is  to  be  expected.  The 
communicative  framework  of  the  design  sup- 
port environment  (for  information  exchange 
both  in  terms  of  input  as  well  as  the  provision 
of  feedback  to  the  user)  should  realistically  con- 
sider the  information-processing  capacity/char- 
acteristics of  the  users  (mainly  practitioners), 
while  maintaining  sufficient  flexibility  within 
the  modeling  environment  to  allow  for  the 
evaluation  of  various  design  options. 

A  review  of  the  literature  on  man-machine 
interface  design  suggests  that  two  vital  princi- 
ples need  to  be  considered,  namely,  consistency 
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and  user-centered  design,  particularly  in  refer- 
ence to  the  importance  of  the  "task  domain" 
and  of  making  the  "mechanisms  of  the  system 
match  the  thought  and  goals  of  the  user" 
(Hutchins  et  al.  1985).  In  the  development  of 
the  META-4  computational  environment,  care- 
ful consideration  was  given  to  the  estab- 
lishment of  an  effective  graphic  interface  for 
structuring  data  input  and  visualization  of  the 
hygro-thermal  process.  Building  performance 
simulation  requires  substantial  input  data  in 
terms  of  material  properties  and  configuration, 
dynamic  external  environmental  conditions 
and  the  desired  internal  environmental  condi- 
tions to  satisfy  functional  requirements  and 
comfort  criteria.  In  META-4,  a  comprehensive 
single-screen  graphic  input  format  has  been  de- 
veloped (cp.  Figure  8  and  Figure  9)  for  both  the 
steady-state  and  the  transient  simulation 
modes.  It  allows  the  user  to  immediately  gain 
an  overview  of  the  design  parameters  and  data 
input  requirements  that  are  supported  by  rele- 
vant built-in  libraries.  Inherent  in  this  ap- 
proach is  the  assurance  of  data  acceptability 
(type  matching)  and  dimensional  accuracy.  The 
graphic  interface  is  designed  to  efficiently  facili- 
tate the  management  and  processing  of  data 
relevant  to  the  application. 

Material  specification  and  the  configuration 
of  any  arbitrary  multi-layered  building  enclo- 


sure component  is  supported  by  a  built-in  mate- 
rial library.  In  the  building  industry,  new  mate- 
rials are  continuously  being  introduced. 
Therefore,  it  is  necessary  for  the  library  of  ma- 
terials to  be  customizable  and  upgradable  (cp. 
Figure  10).  This  spread  sheet  input  format  al- 
lows for  future  inclusion  of  categorical  proper- 
ties for  each  material  that  would  be  required 
for  other  building  performance  simulation, 
such  as  acoustics  and  lighting.  Retrieval  of  pre- 
viously created  construction  files  is  also  facili- 
tated (cp.  Figure  1 1 ). 

The  result  output  interface  tracks  the  "real- 
time" simulation  process  and  presents  the  re- 
sults numerically  and/or  graphically  to  the 
user.  The  benefit  of  the  graphic  output  is  to  en- 
able the  user  to  observe  critical  conditions  that 
may  arise  during  the  simulation  period.  The 
graphical  representation  of  temperature,  water 
vapor  and  moisture  distribution  within  the  en- 
closure component  sub-layers  provides  an  effec- 
tive representation  of  the  pattern  of  behavior  of 
the  component  as  the  boundary  conditions  vary 
in  time  (cp.  Figure  12).  The  ability  to  conduct 
parametric  simulation  studies  and  observe  the 
response  patterns  graphically  promises  to  be 
particularly  effective  from  a  pedagogical  point 
of  view. 
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Figure  8.  Single-screen  data  input  formal'  for  the  steady-slate  (Glaser)  simulation  mode 
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Figure  1  2.  Example  of  a  graphic  output  window  illustrating  saturation  and  "true"  vapor  pressure  curves  as  well  as  condensation  layers 
and  zones. 

(Component  layers:  1 .  Gypsum  Drywall,  2.  Fiberglass,  3.  Aluminum  Foil,  4.  Polyurefhane  Foam  Board,  5.  Building  Paper,  6.  Air 
Space  (non-vent.),  7.  Granite.) 
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CONCLUSION 


The  current  state  of  performance  repre- 
sentation in  computer-aided  architectural  de- 
sign environments  is  characterized  by  a 
dilemma  between  the  two-fold  requirement  of 
rigorous  modeling  capabilities  (crucial  for  sup- 
porting competent  design  decision  making)  and 
integration  (essential  for  supporting  participa- 
tory approaches  to  design  through  effective 
communicative  structures).  An  original  effort 
in  developing  the  META-4  simulation  environ- 
ment has  been  introduced  to  illustrate  the  im- 
portance of  both  rigorous  performance 
representation  and  effective  communicative 
structures.  Preliminary  experimental  applica- 
tion of  the  program  within  the  personal  com- 
puter domain  has  yielded  encouraging  results 
toward  meeting  the  objectives  set  forth  in  this 
research  effort. 
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A  Simple  Algorithm  for  Determining  the  Thermal  Resistance 
of  a  Wall 


Elisabeth  Kossecka,  Jan  Kosny,  and  Jeffrey  E.  Christian 


ABSTRACT 


INTRODUCTION 


This  paper  presents  a  simple  mathematical 
method  of  estimating  wall  thermal  resistance 
based  on  measurements  of  the  wall  surface  tem- 
perature and  heat  flux.  The  article  analyzes  in- 
tegral formulae  for  the  heat  flow  in  finite  time 
across  a  wall  surface  at  prescribed  surface  tem- 
peratures, with  the  terms  related  to  heat  capac- 
ity separated.  The  solution  makes  it  possible  to 
minimize  the  wall  measurement  period  neces- 
sary to  estimate  the  wall  current  time  material 
thermal  resistance  for  time  intervals  with  close- 
to-zero  thermal  flux.  This  method  does  not  de- 
tect dynamic  changes  in  thermal  properties 
such  as  those  caused  by  thermal  drift  or  mois- 
ture adsorption  /  desorption.  The  algorithm 
was  examined  during  field  tests  of  foundation 
wall  systems  conducted  at  Oak  Ridge  National 
Laboratory. 


The  thermal  efficiency  of  passive  solar  sys- 
tems has  been  studied  by  J.  Kosny,  E. 
Kossecka,  and  their  colleagues  as  part  of  the 
Polish  Program  of  Energy  Conservation  in 
Buildings  (Kosny,  Starakiewicz,  and  Szyman- 
ski  1986-90)  (Kossecka  et  al.  1989).  As  part  of 
that  work,  they  have  studied  theoretically  the 
problem  of  heat  exchange  between  a  building 
wall  or  separated  envelope  element  and  its  sur- 
roundings during  long  time  intervals.  This  pa- 
per discusses  an  experimental  application  of  a 
proposed  mathematical  solution  for  the  homoge- 
neous massive  wall  (Kosny  1990)  (Kossecka 
1991)  (Kossecka  1992)  as  a  calculation  method 
increasing  the  accuracy  of  R-value  estimations 
based  on  experimental  results  for  short  meas- 
ure periods  with  the  fluctuating,  close-to-zero 
thermal  flux.  The  accuracy  of  wall  R-value  esti- 
mation during  periodic  changes  in  exterior  tem- 
peratures was  also  discussed  by  Flanders  and 
Mack  (1991).  The  foundation  walls  perform- 
ance investigation  conducted  at  Oak  Ridge  Na- 
tional Laboratory  (ORNL)  by  Jeff  Christian 
helped  verify  the  accuracy  of  this  mathematical 
method  (Christian  1991)  (Christian  and 
Hameister  1993  draft). 
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tional Laboratory  in  Oak  Ridge,  Tenn.,  currently  conduct- 
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HEAT  FLOW  ACROSS  THE  WALL 

Consider  a  plane  wall  of  thickness  L,  in 
which  the  one-dimensional  heat  transfer  condi- 
tion is  satisfied.  Let  Fourier's  law  and  the  one- 
dimensional  heat  equation  be  satisfied 
(Carslaw  and  Jaeger  1959). 


k  dx 


(1) 
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dT      d 

pc¥=^ 


dx 


(2) 


The  thermophysical  properties  of  the  wall — 
thermal  conductivity  k,  specific  heat  c,  and  den- 
sity p — are  constant  in  time.  The  temperature 
in  the  wall  is  represented  by  the  function  T\x,t) 
of  the  spatial  coordinate  x  and  the  time  t.  The 
heat  flux  is  also  represented  by  the  function 
q(x,t). 

In  the  case  of  a  multilayer  wall  composed  of 
materials  of  different  thermal  properties,  Equa- 
tion (2)  has  to  be  understood  symbolically  as  a 
set  of  heat  conduction  equations  written  for  the 
separate  layers  and  the  equations  of  continuity 
of  temperatures  and  heat  flow  rates  at  the  in- 
terfaces. The  following  equations  assume  a  coor- 
dinate system  in  which  wall  surfaces 
correspond  to  the  planes  x  =  0  and  x  =  L 
(Figurel).  The  boundary  values  of  temperatures 
and  heat  fluxes  are  denoted  by  7X0),  7\L),  g(0), 
and  g(L)  respectively;  the  fluxes  in  the  direc- 
tion perpendicular  to  the  wall  surfaces  are  de- 
noted by  qn(Q)  and  qn(L).  In  the  assumed 
coordinate  system: 


Let R0.x  and  RxL  denote  the  heat  resistances 
of  the  wall  layers  enclosed  in  the  intervals 
(0,  x)  and  {x,  L)  respectively: 


R 


0-*- 


dx' 
k(x') 


and 


Rx-L  =  J 


dx' 

kix') 


(4) 


The  total  heat  resistance  is  denoted  by  R. 
R0.x,  Rx.i,  and  R  satisfy  the  following  identity: 
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(>  <  Rx-L 

+  — 7T-  =  1 
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(5) 


After  multiplying  the  conduction  equation 
(Equation  2)  for  temperature  T(x,t)  by  RX.JR 
and  integrating  by  parts,  from  0  to  L,  the  ex- 
pressions appearing  on  its  right  and  left  sides, 
we  can  obtain  the  following  equality: 
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Figure  1 .  Experimental  wall  schematic. 


The  analogous  equations  in  which  q(L)  ap- 
pears are  to  be  obtained  by  multiplying  (Equa- 
tion 2)  by  R0JR: 


<7„(0)  = 


and 
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dT 
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(8) 


Assuming  now  the  constant-in-time  thermo- 
physical properties  of  the  wall  and  integrating 
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overtime  (Equations  7  and  8),  we  obtain  the  fol- 
lowing integral  formulae  for  the  total  amount 
of  heat  flow  across  the  surfaces  of  the  wall: 


Qn  {0,h,t)  =  -^\  df  [7X0)  -T{L)]  - 


i?v_ 


J  dx  -g*  pc[T(x,  0  -  T(x,  t0)] 


R 


and 


(9) 


viously  about  the  time  independence  of  the 
physical  parameters  of  materials,  the  quanti- 
ties Qc0  and  QcL  depend  exclusively  on  the  tem- 
perature differences  in  the  initial  and  final 
states;  if  there  is  no  temperature  difference, 
the  Qc0  and  Qd  vanish.  For  prescribed  tempera- 
tures of  the  wall  surfaces  7X0)  and  T(L),  they 
are  the  only  unknown  components  of  the  heat 
flows  Q„(0)  and  Qn(L)  across  those  surfaces.  If 
there  is  no  temperature  difference  between  the 
surfaces  (symmetric  heating  or  cooling),  Qm 
and  Qci  are  almost  equivalent  to  Q„(0)  and 

Qn(D. 


Qn  (L,  t0  ,  t)  =  -  i  \  df  [T(L)  -7X0)]  - 
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R, 
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(10) 


Denoting  by  T,  T(L),  q  (0),  and  q  (L)  the 
boundary  temperatures  and  heat  fluxes  aver- 
aged over  the  time  interval  of  the  length  At,  the 
formulae  (9)  and  (10)  may  be  rewritten  as  the 
following: 


Summing  up  both  sides  of  Equations  (9)  and 
(10)  and  taking  into  account  identity  (Equation 
5),  we  can  write  the  heat  balance  of  the  wall  in 
the  following  form: 


Qn(0)  +  Qn(L)  =  -Qc 


(11) 


Qn  (0) 


T  (0)  -  T  (L) 
R 


QcO 

At 


and 


Qn  (L) 


T(L)-T  (0)  _  Qcl 
R  At 


(15) 


(16) 


where  Qc  denotes  the  difference  of  the  amount 
of  heat  stored  in  the  wall: 


Qc(t0,t)  =  ldxpc[T(x,t)-T(x,t0)~\   (12) 


Let  Qco,  and  Q<.l  denote  the  following  quanti- 
ties: 

QcQ  =  \dx  -j^pc[T(x,t)-T(x,tQ)l(.m 


Expressions  (15)  and  (16)  can  be  solved  with 
respect  to  VR,  which  gives  the  following  de- 
pendence: 


Qn  (0) 


R 


t [1+Y„(0)]  = 


T(L)  -  T{0) 


^%-[l  +  Y»(L)] 


T(0)  -  T(L) 


where: 


(17) 


and 


R, 


Qd  =  }dx  -jripc[T(x,t)-T{x,t0)]{14) 


R 


Y»(0)  =  ; 


1       Q. 


9n(0) 


At 


and 


(18) 


Qc0  and  Qd  may  be  interpreted  as  the  compo- 
nents of  the  total  heat  flow  across  the  wall  sur- 
faces contributing  to  the  difference  in  the 
amount  of  heat  stored  in  its  volume  of  unit 
cross-section.  Under  the  assumptions  made  pre- 


Y»tfO  =  - 


1      Q, 


Qn(L) 


cL 
At 


(19) 
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HOMOGENOUS  WALL  AND  ASYMPTOTICALLY 
STEADY  HEAT  FLOW 

An  examination  of  Equations  (13)  and  (14) 
indicates  that  the  component  of  Qc  correspond- 
ing to  a  given  surface  is  generally  more  sensi- 
tive to  the  temperature  difference  in  the 
vicinity  of  this  surface  than  in  the  vicinity  of 
the  opposite  one.  However,  there  are  additional 
factors — the  thermal  capacity  and  resistance 
distribution.  In  general,  evaluation  of  the  quan- 
tities Qc0  and  Qd,  contributing  to  the  difference 
in  the  amount  of  heat  stored  in  the  wall  Qc,  for 
prescribed  temperatures  T(0,t)  and  T(L,t), 
makes  it  necessary  to  solve  the  transient  heat 
conduction  problem.  Explicit  expressions  can 
be  written  for  the  heat  conduction  process,  for 
which  the  initial  and  final  temperature  distri- 
bution in  the  wall  correspond  (in  the  sense  of 
asymptotical  convergence)  to  the  steady  state 
heat  flow.  They  can  be  written  as  follows: 


r(x)  =fT(L)  +^r 


Rq-x    _    X_ 

R      ~  L 


T(0) 


(20) 


(21) 


and 


2L_ 


x-L 


L-x 


R 


(22) 


The  amount  of  heat  accumulated  by  the  wall 
may  be  expressed  in  the  following  forms  which 
are  derived  by  integration  Equations  (12),  (13), 
and  (14): 


Qco  =  Lpc 


and 


Qd  =  Lpc  ~ 


|aT(0)+|aT(L) 


|  A  T(0)+|  A  T(L) 


(23) 


(24) 


where 

A  2X0)  =  2X0,  t)  -  2X0,  to)  . 


and 


A  2XL)  =  T(L,  t)  -  T(L,  t0) 


(26) 


(27) 


Qc0  is  twice  as  sensitive  to  the  variation  of 
T(0)  as  to  the  variation  of  T\L);  the  opposite 
holds  for  Qd.  Equations  (23)  and  (24)  may  be 
helpful  in  estimating  the  magnitude  of  the 
terms  QC(jAt  and  Qci/At  in  expressions  (15)  and 
(16)  for  g„(0)  and  qn(L)  and  the  magnitude  of 
the  ratios  y„(0)  and  yB(L)  defined  by  (18)  and 
(19).  This  can  help  in  determining  wall  R-value 
by  using  formula  (17).  When  the  performance 
of  the  massive  storage  wall  with  very  light  or 
transparent  insulation  is  examined,  QcjAt  and 
Qc^/Af  values  may  be  inspected  through  tem- 
perature measurements  with  the  help  of  sen- 
sors located  on  the  interior  and  exterior 
surfaces  of  the  wall,  and  inside  the  wall  (Kosny 
1990). 

EXPERIMENTAL  APPLICATION 

Expression  (17)  may  serve  as  an  algorithm 
for  determining  the  thermal  resistance  of  the 
wall  on  the  basis  of  measurements  of  surface 
heat  fluxes  and  temperature.  For  steady  inside 
surface  temperature  A  T  (L,  t0,  t)  =  0,  yn  (L,  yt) 
can  be  described  by  the  following  dependence: 


yn(L,  At)  = 


(A  2X0) 
At 


A^m-rmi-ASa 


L  pc 


At 


(28) 


and 


where  wall  material  thermal  conductivity  is  de- 
noted by  k. 

Let  A,  B,  C  denote  the  following  quantities: 

A  AT{0) 


Qc  =  Lpc  -  [  A  2X0)  +  A  T(L)]  , 


(25) 


B  =  T(0)  -  T(L)  , 
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6k 
L  pc 


Then  formula  (28)  may  be  rewritten  as 
below: 


(29) 


ln(L,  At)  = 


A 


CB  -  A 


(30) 


In  light  of  the  fact  that  for  long-time  meas- 
urements the  terms  Qc0/At  and  Qa^/At  often 
may  be  neglected,  values  y„(0)  and  y„(L)  are 
very  small  compared  with  unity.  That  is  why 
approximate  wall  thermal  resistance  R°,  in 
case  of  simple  averaging  calculations,  can  be  es- 
timated from  the  following  expression: 


_1_ 

R 


Qn(0) 


Qn(D 


T  (L)  -  T(0)       T  (0)  -  T(L) 


(31) 


If  Equations  (16)  and  (31)  are  compared,  the 
following  dependence  can  be  obtained: 


^  =  1  +  yn(JL,At) 


(32) 


Equation  (31)  is  commonly  used  to  estimate 
approximate  wall  R-values  based  on  surface 
thermal  flux  and  temperature  measurements. 
According  to  Equation  (32),  the  error  of  such  es- 
timation is  as  small  as  the  yn(L)  value  which  is 
close  to  zero.  That  is  why  a  knowledge  of  the 
Y„(L)  value  (referred  to  hereafter  as  a  correc- 
tion ratio)  may  help  to  estimate  an  accuracy  of, 
based  on  Equation  (30),  approximate  R-value 
calculations.  In  Figure  2,  values  of  the  correc- 
tion ratio  are  shown  for  several  solid  concrete 
walls,  12,  8,  and  4  inches  thick,  made  of  120 
lb/ft3  and  40  lb/ft3  concretes.  It  can  be  seen  that 
for  surface  temperature  differences  T(Q)-T(L) 
equal  50°F,  a  value  of  yn(L)  can  be  smaller  than 
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Figure  2.  Correction  ratios  for  solid  concrete  walls. 
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0.05  for  4-day  measure  periods  in  case  of  12- 
inch  walls  and  for  2-day  periods  for  8-inch 
walls.  In  case  of  4-inch  walls,  a  value  of  7„(L) 
can  be  smaller  than  0.01  for  2-day  measure  pe- 
riods, and  for  the  same  T(0)-T(L),  surface  tem- 
perature differences.  A  reason  for  achieving 
such  small  values  of  yn(L)  lies  in  the  fact  that 
vary  small  values  of  differences  were  assumed 
between  initial  and  final  wall  surface  tempera- 
tures, 7X0). 

In  practice,  based  on  temperature  profile  ob- 
servations, it  is  possible  to  select  time  periods 
where  the  differences  between  initial  and  final 
temperatures,  A7X0),  are  minimal,  and  thermal 
flux  does  not  change  in  direction  of  flow.  Figure 
3  presents  values  of  temperature  differences 
A7X0)  that  can  lead  to  estimation  R-values  with 
an  error  of  less  than  5  percent  (correction  ratio 
values  yn  are  less  than  0.05). 

Based  on  the  ORNL  foundation  wall  experi- 
mental results  (Christian  and  Hameister  1993 
draft),  the  authors  selected,  for  an  uninsulated 
wall,  measurement  periods  where  the  initial 
and  final  states  of  the  wall  were  almost  equal 
to  minimize  A7X0).  Then  R-values  were  calcu- 


lated. By  selecting  time  intervals  during  which 
the  initial  and  final  temperatures  of  the  wall 
surface  where  the  thermal  flux  meter  was  lo- 
cated, were  almost  equal,  it  was  possible  to  esti- 
mate R-values  even  for  2-day  measurement 
periods  (Figure  4).  Results  were  compared  with 
R-values  obtained  by  use  of  the  PROPOR  com- 
puter code  (Beck,  Petrie,  and  Courville  1991) 
and  with  R-values  calculated  according  to  Equa- 
tion (31),  both  for  7-day  time  intervals.  These 
are  shown  in  Figure  5.  Most  of  the  results  for 
the  spring  and  fall  months  are  in  good  agree- 
ment. For  summer,  because  of  very  small  and 
fluctuating  thermal  flux,  calculated  wall  R-val- 
ues are  divergent.  It  was  possible,  by  selecting 
time  intervals,  to  obtain  relatively  accurate 
wall  R-values  even  for  2-day  measurement  peri- 
ods (too  short  for  the  traditional,  averaging 
way  of  R-value  computations,  and  for  PROPOR 
analyses).  For  selected  time  periods,  it  was  also 
possible  to  estimate  wall  R-values  for  consider- 
ably smaller  values  of  average  thermal  flux.  Ta- 
ble 1  presents  calculated  wall  R-values  for  June 
1991  obtained  by  using  the  PROPOR  computer 
code  and  Equation  (31)  for  selected  shorter  peri- 
ods. 
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time  [h] 


Figure  3.  Maximum  values  of  AT(0)  for  concrete  solid  walls  giving  correction  ratios  less  than  0.05. 
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CONCLUSIONS 


From  the  analysis  of  the  formulas  (9),  (10), 
(15),  (16),  (17),  (18),  and  (20),  a  series  of  conclu- 
sions may  be  developed  that  can  be  useful  in  ex- 
perimental analyses  of  wall  thermal  properties. 
Some  of  them  are  as  follows: 

•     The  heat  flow  through  the  wall  is 
the  same  for  all  heat  transfer  proc- 
esses with  the  same  initial  and  final 
surface  temperatures,  and  thermal 
flux  not  changing  in  direction  of 
flow. 


Over  long  time  intervals,  the  mass 
distribution  inside  the  wall  has  no 
influence  on  the  heat  flow  under  cy- 
clical surface  temperatures. 

When  the  thermal  resistance  of  the 
wall  is  to  be  determined  on  the  ba- 
sis of  surface  temperatures  and 
heat  flux  measurements,  with  the 
use  of  approximate  formula  (19), 
the  error  can  be  minimized  if  the  ex- 
periment is  conducted  so  that  the  in- 
itial and  final  temperatures  at  the 
point  of  the  thermal  flux  measure- 
ment are  close.  This  can  shorten 
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Figure  4.  R-value  calculations  for  periods  selected  by  thermal  flux  and  temperature  profile  observations. 
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measurement  periods,  which  aids 
the  experimental  estimation  of  wall 
thermal  properties  for  time  inter- 
vals with  small  and  fluctuating  ther- 
mal flux. 

The  integral  formulae  (9)  and  (10) 
describe  the  total  amount  of  flow 
across  the  surfaces  of  the  wall.  They 
can  be  used  in  analyzing  the  influ- 
ence of  variations  of  the  tempera- 
ture profile  inside  the  wall  on  the 
heat  flow  across  its  surfaces  for  a 
given  wall  resistance  and  thermal 
capacity  distribution. 


•     The  correction  factor  Yn,  describing 
inaccuracy  of  R- value  estimation  by 
approximate  Equation(31),  may  be 
predicted  by  use  of  roughly  esti- 
mated wall  parameters  and  the  for- 
mulae for  the  steady  state  heat  flow. 

Some  of  the  above  conclusions  are  intuitively 
obvious  or  well  known;  however,  it  is  worthy  to 
note  that  this  paper  presents  a  strict  mathe- 
matical basis  for  them  in  the  form  of  the  inte- 
gral formulae  for  heat  flow  across  wall 
surfaces. 
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Figure  5.  Experimental  R-values  for  several  calculation  methods. 
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Testing 
interval 


Table  1 .  Comparison  of  calculated  R-values 


Length  of 
interval  (h) 


Average 
thermal  flux 
[BtuUvft2)-11 


Calculated 
R-values 

[h-ft2-F/Btu] 


Using 
PROPOR 


Using 
Eq.31 


Entire  period 
Hours  51-96 
Hours  35-87 


Period  beginning  June  3,  1991 

96  1.029375 


45 
52 


1.332706 
1.427924 


0.795548  0.862072 
0.830663 
0.789244 


Entire  period 
Hours  1-76 
Hours  14-60 


Period  beginning  June  11,  1991 

96  0.738645 


75 
46 


0.855394 
0.855744 


1.109011 
1.033225 
0.992541 


Entire  period 
Hours  2-54 
Hours  18-61 


Period  beginning  June  17,  1991 

144  0.382013 


52 
43 


0.523018 
0.495454 


1.271223 
1.096681 
1.171559 


Entire  period 
Hours  16-116 
Hours  5-122 


Period  beginning  June  24,  1991 

168  0.812321 


100 

117 


1.026039 
0.948389 


0.936689 
0.861237 
0.877580 
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Geosynthetics  in  Residential  Construction 

John  S.  Horvath,  Ph.D.,  P.E. 
Associate  Professor  of  Civil  Engineering 
Manhattan  College 
Bronx,  New  York 

The  civil  engineering  profession  is  starting 
to  take  greater  interest  in  applying  geosynthet- 
ics technology  to  residential  and  light  commer- 
cial/industrial structures.  The  general  neglect 
of  this  market  to  date  is  based  primarily  on  the 
educational  emphasis  to  which  all  civil  engi- 
neers are  exposed.  By  omission  from  course  ma- 
terial, civil  engineers  develop  the  perspective 
that  residential  and  similar  construction  does 
not  require  engineering  attention.  As  the  read- 
ers of  the  Building  Research  Journal  are  well 
aware,  this  is  a  fallacy  that  needs  correction. 
The  purpose  of  this  article  is  to  offer  an  over- 
view of  some  developments  that  the  writer  is  fa- 
miliar with  that  have  particular  application  in 
residential  and  related  construction.  It  is  hoped 
that  this  will  initiate  a  discussion  among  the  di- 
verse interests  who  are  involved  in  this  con- 
struction area. 

A  concept  that  is  currently  in  the  early 
stages  of  development  is  that  of  a  "geosynthetic 
envelope"  around  the  below-ground  portion  of 
the  structure.  Among  the  various  geosynthetic 
products  that  are  available,  the  two  that  are 
the  basic  elements  in  this  concept  are  geofoams 
and  drainage  geocomposites.  The  concept  is  to 
place  layers  of  these  two  materials  around  the 
entire  below-ground  portion  of  a  structure  with 
the  drainage  geocomposite  being  the  outer 
layer  that  is  in  contact  with  the  ground. 

Geofoam  is  a  new  but  accepted  term  for  any 
foam  product  used  in  a  below-ground  applica- 
tion. Based  on  experience  to  date,  the  primary 
geofoam  materials  are  the  two  types  of  rigid  cel- 
lular polystyrene  (RCPS):  expanded  polysty- 
rene (EPS  or  "headboard")  and  extruded 
polystyrene  (XPS).  Conventionally,  RCPS  is  in- 
stalled as  panels  typically  in  the  range  of 
25mm  to  75mm  thick  placed  on  the  exterior  of 
the  below-ground  walls  and,  in  some  cases,  be- 
neath the  basement  slab.  However,  there  has 
been  a  rapid  growth  in  systems  where  RCPS 
panels  or  hollow  blocks  are  used  as  formwork 


for  a  poured-in-place  concrete  wall.  One  of  the 
growing  concerns  with  using  RCPS  in  below- 
ground  applications  is  damage  caused  by  insect 
infestation.  A  relatively  recent  development  is 
a  passive  treatment  that  is  applied  during 
manufacture.  Research  indicates  that  this 
treatment  is  successful  in  deterring  insect  infes- 
tation, at  least  for  termites  and  carpenter  ants. 
At  the  present  time  this  treatment  is  available 
only  in  North  America  and  only  for  EPS  panels 
from  certain  manufacturers. 

Traditionally,  the  geofoam  layer  is  envisaged 
only  as  providing  the  function  of  thermal  insu- 
lation to  reduce  life-cycle  energy  costs.  In  this 
regard  it  is  worth  noting  that  not  only  can  ther- 
mal insulation  reduce  costs  during  the  heating 
season  but  there  is  also  a  potential  for  reducing 
dehumidification  costs,  at  least  for  below- 
ground  space,  during  the  summer.  This  is  be- 
cause the  presence  of  thermal  insulation  moves 
the  dew  point  from  inside  the  below-ground 
space  to  outside.  Electric-energy  cost  savings 
from  reduced  dehumidification  should  there- 
fore be  considered  where  appropriate.  A  new 
thermal  insulation  application  in  the  U.S.  is 
called  the  "frost-protected  shallow  foundation" 
(FPSF)  concept.  This  applies  to  structures  in 
relatively  cold  climates  that  do  not  have  below- 
ground  space.  In  these  areas,  the  structure  foot- 
ings must  be  founded  at  significant  depths, 
sometimes  in  excess  of  two  meters,  to  be  below 
the  seasonal  frost  depth.  This  adds  a  signifi- 
cant cost  to  construction.  By  insulating  the 
structure  perimeter,  footing  embedment  can  be 
reduced  to  as  little  as  500mm.  In  addition, 
there  will  be  life-cycle  energy  savings  as  well. 
This  concept  has  been  used  in  Europe  and  else- 
where for  more  than  30  years.  Current  U.S.  ac- 
tivities involve  some  instrumented  houses  in 
the  northern  U.S.  as  well  as  efforts  to  incorpo- 
rate provisions  for  such  construction  in  the  vari- 
ous model  building  codes  in  use  throughout  the 
U.S. 
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Some  new,  additional  functions  of  geofoam 
are  currently  being  investigated.  This  involves 
using  the  foam  as  a  "compressible  inclusion"  to 
reduce  earth  pressures  acting  on  below-ground 
walls,  especially  in  situations  involving  seismic 
loading,  swelling  clays,  and  freezing  soils.  The 
potential  benefit  of  these  new  applications  is 
significant  especially  considering  that  the  geo- 
foam will  still  provide  traditional  value  as  ther- 
mal insulation.  At  the  present  time,  research 
indicates  that  only  panels  of  EPS  can  be  manu- 
factured to  provide  the  necessary  properties  to 
act  as  a  compressible  inclusion.  This  is  fortui- 
tous in  view  of  the  fact  that  only  EPS  panels 
can  be  treated  against  insect  infestation  as 
noted  above. 

The  use  of  a  drainage  geocomposite  can  also 
provide  several  benefits.  In  addition  to  the  obvi- 
ous one  of  intercepting  ground  water  before  it 
enters  below-ground  space,  a  drainage  layer 
minimizes  the  water  absorbed  by  the  geofoam 
layer  thereby  maximizing  its  thermal  effi- 
ciency. There  is  also  a  potential  that  the  drain- 
age layer  can  be  used  as  part  of  a  system  to 
intercept  ground-borne  gases  such  as  radon. 
The  use  of  a  geocomposite  as  the  drainage  ma- 
terial offers  benefits  over  the  use  of  natural  ag- 
gregates (gravel  or  crushed  rock)  both  in  terms 
of  cost  in  many  areas  as  well  as  construction 
cost  savings  in  ease  of  placement.  There  is  even 
one  drainage  geocomposite  that  has  a  signifi- 
cant thermal-insulative  value  that  can  contrib- 
ute in  this  area  as  well. 

In  summary,  the  use  of  a  "geosynthetic  enve- 
lope" around  the  below-ground  portion  of  struc- 
tures has  great  potential  for  providing 
simultaneous,  multiple  benefits  in  a  cost-effec- 
tive manner.  I  welcome  contact  and  comments 
from  those  with  similar  research  interests.  I 
can  be  reached  at: 

Prof.  John  S.  Horvath,  P.E. 

Associate  Professor  of  Civil  Engineering 

Manhattan  College 

Civil  Engineering  Department 

Bronx,  New  York  10471 

Telephone:  718-920-0177 

Fax:  718-796-9812 


From  Other  Universities  — 


University  of  Illinois  at  Urbana-Champaign 
School  of  Architecture-Building  Research  Council 
(BRC) 

The  Council  is  developing  its  role  as  a  na- 
tional publisher  of  research-  and  application- 
oriented  materials.  Projects  include  publication 
of  residential  home-owner  materials;  transla- 
tion of  selected  materials  into  Spanish,  and  a 
national  solicitation  of  manuscripts  appropri- 
ate to  research,  technology,  and  home-owner 
publication. 

BRC  research  efforts  include  an  agreement 
with  State  Farm  to  investigate  the  problem  of 
burst  water  pipes  resulting  from  subfreezing 
temperatures.  CertainTeed  will  continue  and 
expand  its  support  of  attic  performance  re- 
search. The  project,  "Design  and  Implementa- 
tion of  a  Home  Energy  Rating  System", 
includes  information  transfer  to  the  home- 
building/sales  industry;  consumer  education; 
home  energy  rater  training  and  support;  and  in- 
formational data-basing. 

Through  seminars  sponsored  by  the  Illinois 
Attorney  General's  office,  BRC  provided  techni- 
cal information  for  the  repair  of  flood-damaged 
houses  to  victims  of  the  1993  Mississippi  flood. 
These  seminars  provided  citizens  with  con- 
sumer protection  information.  Other  partici- 
pants included  the  Illinois  chapter  of  NAHB 
and  FEMA. 

Due  to  staff  cutbacks,  increased  research  op- 
portunities, and  spiraling  demands  for  public 
service  during  the  1993  summer  flood,  BRC 
was  forced  to  suspend  its  housing  advisory  line 
service.  Architects,  contractors,  and  home  own- 
ers in  Illinois  and  across  the  nation  have  come 
to  depend  on  this  service  to  answer  home  con- 
struction, design,  maintenance,  and  remodeling 
questions.  Funds  are  being  sought  to  hire  a 
housing  expert  half-time  to  reinstate  the  line. 

1994  is  the  Council's  fiftieth  anniversary  and 
the  Year  of  the  Alumni.  BRC  encourages  all 
who  have  worked  with  researchers  at  the  Coun- 
cil to  contact  us  by  calling  us  at  800-336-0616, 
or  by  writing  to  us  at  BRC  / 1  East  St.  Mary's 
Road/  Champaign,  IL  61820.  We  would  like  to 
hear  where  you  are  and  what  you  are  doing, 
whether  you  would  like  a  newsletter  about 
Council  and  alumni  activities,  and  how  you 
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would  like  to  renew  ties  to  BRC  and  its  mis- 
sions. 

University  of  Florida 

Shimberg  Center  for  Affordable  Housing 

The  Shimberg  Center  organized  the  4th 
Rinker  International  Conference  on  Building 
Construction  with  the  theme  of  "Affordable 
Housing:  Present  and  Future."  As  a  result  of 
this  conference,  the  Center  established  an  inter- 
national network  of  over  100  housing  and  hous- 
ing-related researchers  and  professionals  in 
over  50  countries,  plus  a  large  number  in  the 
U.S.  This  network,  in  turn,  led  the  Center  to 
propose  re-establishment  of  W63:  Low-cost 
Housing,  one  of  the  Working  Commissions  of 
the  Netherlands-based  International  Council 
for  Building  Research  Studies  and  Documenta- 
tion (CIB).  The  success  of  the  conference  and 
the  subsequent  national  and  international  rec- 
ognition was  directly  related  to  the  support  re- 
ceived by  the  Center  from  Fannie  Mae's  Office 
of  Housing  Research,  the  National  Housing  En- 
dowment, and  the  Nations  Bank. 

University  of  Oregon 

The  Center  for  Housing  Innovation  (CHI) 

CHI  has  undertaken  the  Energy  Efficient  In- 
dustrialized Housing  Research  Program 
(EEIH).  This  program,  funded  by  the  U.S.  De- 
partment of  Energy,  private  industry,  state  gov- 
ernments and  utilities,  addresses  the  need  to 
increase  the  energy  efficiency  of  industrialized 
housing  through  both  design  and  production  im- 
provements. The  Center's  goal  is  to  conduct  re- 
search that  will  develop  techiques  to  produce 
marketable  industrialized  housing  that  is  25 
percent  more  energy  efficient  than  the  most 
stringent  U.S.  residential  codes  now  require 
and  will  cost  less. 

University  of  Georgia 
Housing  Research  Center 

The  Center  completed  work  on  "An  Assess- 
ment of  Changes  in  Energy  Behavior:  Revis- 
ited", funded  by  the  Governor's  Office  of 
Energy  Resources,  and  "An  Analysis  of  Water 
Quality  from  Private  Wells  in  Georgia:  The  In- 
fluences of  Land  Use,  Enonomic  and  Demo- 
graphic Factors  and  Strategies  for  Families  at 
Risk",  which  was  funded  by  the  USDA  Agricul- 
tural Research  Service. 


Texas  A&M  University 

Center  for  Housing  Research  and 

Urban  Development 

The  Center  completed  two  community  re- 
source centers  in  the  colonias  along  the  Texas- 
Mexico  border.  Centers  in  the  Cameron  Park 
colonia  near  Brownsville,  and  a  colonia  in  Pro- 
greso,  near  Weslaco,  Texas,  were  completed 
during  December,  1993.  Ribbon  cutting  ceremo- 
nies are  being  planned  for  the  Laredo  and  El 
Paso  areas  and  sites  for  the  fifth  and  sixth  cen- 
ters are  being  investigated. 

Colorado  State  University 
Housing  Research  Center 
Program  in  Computer-Integrated 
Construction  Management 

The  current  focus  areas  in  the  Center  and 
the  Program  include  corporate  and  project  man- 
agement, computer-integrated  design  and  con- 
struction planning,  international  construction 
management,  stabilization  and  historic  preser- 
vation, and  indoor  air  quality. 

The  Center  and  the  Program  have  been  in- 
volved in  the  following  initiatives  and  projects 
related  to  international  construction  manage- 
ment: 

•  Construction  safety  training  for  the 
Ministry  of  Labour,  Cyprus; 

•  Computer-integrated  planning  and 
scheduling  training  for  officials 
from  Pakistan; 

•  Provision  of  environmental  techni- 
cal assistance  as  part  of  a  team  ef- 
fort to  Ex-Soviet  Union  countries; 

•  A  visit  by  the  former  president  of 
the  School  of  Planning  and  Architec- 
ture, New  Delhi,  India,  to  explore 
the  possibilities  of  mutual  coopera- 
tion. (The  School  has  graduate  pro- 
grams in  Building  Engineering  and 
Management  and  Housing  Planning 
Development.) 

•  A  seminar  presentation  on  "Com- 
puter-Integrated Design  and  Con- 
struction" at  the  School  of  Planning 
and  Architecture,  New  Delhi,  India. 
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